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Abstract 



At strong magnetic fields double-layer two-dimensional-electron-gas sys- 
tems can form an unusual broken symmetry state with spontaneous inter- 
layer phase coherence. In this paper we explore the rich variety of quantum 
and finite-temperature phase transitions associated with this broken symme- 
C^ I try. We describe the system using a pseudospin language in which the layer 

degree-of-freedom is mapped to a fictional spin 1/2 degree-of-freedom. With 
this mapping the spontaneous symmetry breaking is equivalent to that of a 
spin 1/2 easy-plane ferromagnet. In this language spin-textures can carry 
a charge. In particular, vortices carry ibe/2 electrical charge and vortex- 
antivortex pairs can be neutral or carry charge ite. We derive an effective 
low-energy action and use it to discuss the charged and collective neutral 
excitations of the system. We have obtained the parameters of the Landau- 
Ginzburg functional from first-principles estimates and from finite-size exact 
diagonalization studies. We use these results to estimate the dependence of 
the critical temperature for the Kosterlitz-Thouless phase transition on layer 
separation. 
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I. INTRODUCTION 

Technological progress has made it possible to produce double-layer two dimensional elec- 
tron gas systems of extremely high mobility. As illustrated schematically in Fig.(|I|), these 
systems consist of a pair of 2D electron gases separated by a distance d so small {d ~ lOOA) 
as to be comparable to the typical spacing between electrons in the same layer. In a large 
magnetic field, strong correlations between the layers have long been expected to lead to 
novel fractional quantum Hall effects. Correlations are especially important in the strong 
magnetic field regime because all electrons can be accommodated within the lowest Landau 
level and execute cyclotron orbits with a common kinetic energy. The fractional quantum 
Hall effect occurs when the system has a gap for making charged excitations, i.e. when the 
system is incompressible, and theory has predictedou that at some Landau level filling fac- 
tors, gaps occur in double-layer systems only if interlayer interactions are sufficiently strong. 
These theoretical predictions have recently been confirmeda. More recently work from sev- 
eral different points of viewUB has suggested that inter-layer correlations can also lead to 
unusual broken symmetry states with spontaneous phase coherence between layers which 
are isolated except for inter-layer Coulomb interactions. We have recently arguedEl that it 
is spontaneous interlayer phase coherence which is responsible for the recently discoveredEJ 
extreme sensitivity of the fractional quantum Hall effect at total Landau level filling factor 
z/ = 1 to small tilts of the magnetic field away from the normal to the layers, {u = N/N^ 
where N is the number of electrons and N^ is the number of single-particle levels per Landau 
level.) 

We have previouslyiy presented a phenomenological theory of the rich zero-temperature 
phase diagram associated with spontaneous interlayer coherence. In the present paper we 
provide a detailed microscopic derivation of the effective action used in our phenomenolog- 
ical theory. We also discuss the low-energy neutral and charged excitations of the system 
in some detail. Much of our discussion of the properties of double-layer systems with spon- 
taneous interlayer coherence will be couched in language based on a simple mappingQ of 
the layer degree-of-freedom in a double-layer system to an artificial 'pseudospin' degree-of- 
freedom. In this language the spontaneous-interlayer-coherence broken symmetry appears 
as easy-plane ferromagnetism. The mapping is convenient because the Hamiltonian of the 
system may be simply expressed in terms of pseudospin operators and because some as- 
pects of the physics are familiar when expressed in this way. This mapping is discussed in 
detail in Section ||. In Section [H we present a microscopic derivation of the connection 



between spin-textures and Coulomb charges in our model. This relationship was discussed 
previously by Sondhi et a/.t3 in the context of a Chern-Simons effective field theory for the 
case of a single- layer system at z/ = 1 with weak Zeeman coupling (to the real spin). In 



Section |^ we derive an effective action which describes the low-energy physics of the system 
whenever the system has spontaneous interlayer coherence and is incompressible. We use 
the spin-charge connection and our effective action in Section to discuss the low-lying 
excitations of the system which are formed from vortices in the pseudospin configuration. 
We show that vortices carry charge ±e/2 and that the vortices appear in four fiavors cor- 
responding to the independently available sign choices for vorticity and Coulomb charge. 
Neutral excitations of finite energy can be formed from vortex pairs with both opposite 
vorticity and opposite charge. Collective spin-wave-like excitations of the system also occur 



and dominate response-functions at long-wavelengths. These collective modes are discussed 



in Section 0. We believe that a double-layer system with spontaneous interlayer coherence 
will have a finite-temperature Kosterlitz-Thouless phase transitionO'tl which we discuss in 



Section |VI1] . In the low-temperature phase a kind of superconducting behavior will occur 
in which the linear resistivity vanishes when opposite currents are carried in the two layers. 
Fully microscopic calculations for the double-layer systems, using exact diagonahzation stud- 



ies of finite-size systems and many-body perturbation theory, are reported in Section |V111 



These calculations allow us to estimate the parameters of the low-energy effective action 
and hence to provide quantitative estimates of the dependence of the Kosterlitz-Thouless 
temperature on the separation between the layers. Section |I^ gives a brief summary of the 
Chern-Simons effective field theory description of double layer systems. Finally in Section 
^ we briefiy summarize our findings. A companion paperEj will discuss issues which arise 
when a weak symmetry-breaking tunneling term is added to the Hamiltonian of the double- 
layer sy_stem, particularly those issues which arise from the recent experiments of Murphy 
et. al 



II. SPIN ANALOGY 



I 



We wish to show that the double-layer system at certain total filling factors, particularly 
at u = 1, can be viewed as an easy-plane quantum itinerant ferromagnet. In this section 
we will give a qualitative introduction to the essential ideas of the physical picture. The 
mathematical details of the microscopic physics will be presented in the subsequent sections. 

We will use a 'pseudospin' magnetic language in which pseudospin 'up' ('down') refers 
to an electron in the 'upper' ('lower') layer.Gil Using this language and building upon recent 
progress in understanding the case of single-layer systems at z/ = 1 with real spinE3ll3 we will 
explore the consequences of the mixing of charge and pseudospin degrees of freedom and 
discuss the rich variety of phase transitions controlled by temperature, layer separation, layer 
charge imbalance, and in the companion papert3 we explore the effects of tunneling between 
layers, and magnetic field tilt angles. The present section will be devoted to development 
of a physical picture of the rather counterintuitive concept of spontaneous phase coherence 
between the layers, which in the magnetic analogy corresponds to spontaneous pseudospin 
magnetization. Technical details of the microscopic calculations on which this picture is 
based will be presented in the subsequent sections. 

It is helpful to begin study of the pseudospin analogy by reminding ourselves of the 
unusual properties of a single layer system at i/ = 1 in the limit of zero Zeeman splitting 
(for the 'real' spins) .Ej'Ej In the presence of Coulomb repulsion between the particles, Hund's 
rule would suggest that the system could lower its interaction energy by maximizing its total 
spin since states with maximum total spin are symmetric under spin exchange and hence 
the spatial wave function is necessarily fully antisymmetric. In an ordinary ferromagnet the 
Hund's rule tendency to maximize the total spin is partially counteracted by the increase 
in kinetic energy (due to the Pauli principle) that accompanies spin polarization. In the 
lowest Landau level however, the kinetic energy has been quenched by the magnetic field 
and the system will spontaneously develop 100% polarization. An explicit microscopic wave 
function believedO to exactly describe the ground state of N electrons at z/ = 1 isla 



^ = ^v-| TTTTTTTTTTTT ••• T), (i) 

where ^y is a Vandermonde determinant wave functiontl of the form 

^y = n {z, - Zj)X{eM-\zk?/^^^). (2) 

i<j k 

where (. = {hc/eBY^'^ is the quantized cyclotron orbit radius of the lowest Landau level. The 
first term in Eq.(|ip is simply the Laughlin spatial wave function for the filled Landau level 
and the second term indicates that every spin is up. This state has total spin S = N/2 and 
5*^^ = (A^/2)\l'. Because Coulomb interactions do not directly affect the spins (magnetism 
is caused by Coulomb forces plus the Pauli principle, not magnetic forces!), [H, S^^] = and 
^ is simply one of a total of 25" + 1 degenerate states, all with S = N/2. The other states 
are simply created using the total spin lowering operator S~ = Y^jLi sj which is itself fully 
symmetric under spin exchange. Since the exact ground state (at S^ = S = N/2) is a single 
Slater (Vandermonde) determinant the two-particle distribution function for these states is 
readily computed 

^(|r-r'|) = l-e-'-r'l^/^ (3) 

Here and in the rest of the paper we set the magnetic length, £, to unity. One clearly sees, 
in this expression, the 'exchange hole' which surrounds each particle and lowers its Coulomb 
energy by an amount 

E^ = ^j£r^- \g{r) - 1] = e^^ ~ 64K ^B[Tl (4) 

where the dielectric function of the semiconductor environment is implicit, B\T] indicates 
the magnetic field in Tesla, and the numerical estimate in Kelvins is for the case of GaAs. 
For excited states, the spin wave function is not fully symmetric and thus the spatial wave 
function is not fully antisymmetric. There is then a finite amplitude for particles of opposite 
spin to approach each other closely. The low- lying excited states are 'magnons' and just as 
for the ferromagnetic Heisenberg model on a lattice, the exact single magnon excitations for 
this (itinerant.) magnet can be found. They are labeled by a conserved momentum k and 
have the formEj 

*k = ^^, (5) 

where 



N 



S^^Y^e^-^sJ, (6) 

i=i 

is the Fourier transform of the local spin lowering operator and the overbar indicates pro- 
jection onto the spatial wave functions of the lowest Landau level. The dispersion of these 
excitations is similar (at long wavelengths) to those of the Heisenberg model on a lattice 
and is given byli3£2l 

u^ik) = I -0^V{q) exp(-|g|V2) [1 - cos(z ■ q x k)] , (7) 




where V{q) is the Fourier transform of the electron-electron interaction. For large wave 
vectors one can show that this excitation crosses over from being a collective spin wave 
mode to a single-particle type excitation consisting of a magnetic 'exciton' which is a bound 
state of a spin-flipped particle-hole pair.By 

In the discussion above, we have focussed on the ground state which has definite total 
S and have examined a basis in which 5*^ is a good quantum number. It is convenient for 
later use however to recall that since the Hamiltonian is invariant under spin rotations, we 
could have chosen the spin quantization axis along any direction. Suppose for example we 
considered the ferromagnetic state |\I^(^) in which every electron has the spinor 

(8) 

This state is an exact ground state of the Hamiltonian but has indefinite S^; that is, it is 
made up of a linear combination of all the degenerate S^ eigenstates. The mean value of S^ 
is zero 

{^^\S^\^^) = 0, (9) 

but, 

(^^|S|^^) = - [cos((/^)x + sm{ip)y] , (10) 

showing that the spin is fully aligned but now lies in the xy plane. Fluctuations in S^, while 
non-zero, are relatively small (in the limit of large A^) 

{^.\[S'f\%f" = lNy'. (11) 

This state can be represented as a coherent superposition of the eigenstates of S^ (obeying 
S^lm) = mlm)). For large A^, we have, to a good approximation 

N/2 

1^^) = J2 e-^"''e-'""^\m). (12) 

m=-N/2 

Such a state has a coherence which is analogous to that in the BCS state of a superconduc- 
tor (with S^ playing the role of number operator). There is a definite phase relationship 
between states with different values of S^. Notice that this is not a direct consequence of 
the dynamics, but merely a result of our choice of linear combination of the degenerate basis 
vectors. That is, there are no terms in the Coulomb Hamiltonian which flip spins and yet in 
the ground state there can be a definite phase relationship between amplitudes for different 
numbers of flipped spins. 

We turn now from the application of these ideas to single-layer systems with 'real' spin to 
the analogous ideas for double-layer systems described by pseudospin. We will ignore 'real' 
spin, assuming it to be frozen out by the Zeeman energy, although this is not necessarily a 
valid assumption at low B fields. The spinors 




(13) 

describe states in which the electron is in the upper or lower layer respectively. Thus the 
layer number difference operator is simply the z component of the total pseudospin 

N^-N^ = 2S\ (14) 

One's first reaction upon thinking about pseudospin is that it is a perfectly sensible concept 
as long as it is an Ising-like variable; i.e., each electron is either in the upper layer or in 
the lower, but not both. However it is a fundamental feature of quantum mechanics that 
it is perfectly sensible to talk about states in which there is a coherent superposition of 
two amplitudes and the layer index is therefore uncertain. For certain filling factorsH^S 
the ground state spontaneously develops interlayer coherence and the electron pseudospin 
condenses into a state 

(15) 

magnetized in the xy plane. Such a state has {N^ ~ ^i) = which reduces the charging 
energy of the double layer system. Such a state also has good exchange energy because, if 
two electrons of the same pseudospin orientation (phase (p) approach each other, the spatial 
part of the wave function must vanish. It is this exchange effect which gives rise to the 
finite pseudospin 'stiffness' which polarizes the pseudospins. Of course, in the absence of 
tunneling, Q^ = A^^i — A^^j^ is a good quantum number while our variational wave function 
\1/^ has N^^'^ fluctuations in Q^. In analogy with coherent BCS states however, this is not 
important to the physics (usually), it is simply mathematically convenient not to project 
\I'^ onto a state of definite Q^. Of course it is perfectly possible to do so using 




^Q_=Jdipe-''^-^^^, (16) 

without seriously modifying the good correlations built into the wave function. 

Even though the total number of electrons in a given layer is a good quantum number, 
the dynamics will enforce a definite phase relationship among states with different 5*^ (layer 
charge difference) due to a spontaneously broken U(l) symmetryoQ'S corresponding to ro- 
tations in pseudospin space about the z axis. This is quite analogous to what happens in 
superconductors. For finite layer separation, the charging energy will limit the fluctuations 
in Q- and modified correlations will have to be built into the wave function as discussed in 
Section 0. 

The above discussion for z/ = 1 is, of course, somewhat over-simplified for general filling 
factors. Hund's rule, which suggests that the ground state should have the maximum total 
spin quantum number consistent with the Pauli exclusion principle, does not always apply 
to two-dimensional electrons in the strong magnetic field limit. In particular it is known 
both from theoretical worka'Ejcll and from experimental workB3 that at some filling factors 
incompressible ground states can occur which are spin singlets. c3 We focus our discussion 
here on the case u = 1 where the consequences of spontaneous-interlayer-coherence are 
likely to be most easily observable. At this filling factor there is ample evidence that Hund's 



rule does give the correct answer for the spin quantum number of the ground state. We 
beheve that much of the physics we discuss will occur at any filling factor for which the 
corresponding spin system has an incompressible ground state which is not a spin-singlet. 
In particular, the ground state for u = 1/m for any odd integer m is believed to have 
S = N/2 for any physically realistic interaction. The ground state orbital wavefunctions at 
these filling factors are well approximated by the simple Jastrow wavefunctions discovered 
by LaughlinlHaEa and we will present some results for m = 3 and m = 5 based on these trial 
wavefunctions. 



III. SPIN-CHARGE RELATION 

A. Review of projection onto the lowest Landau level 

A convenient formulation of quantum mechanics within the subspace of the lowest Lan- 
dau level (LLL) was developed by Girvin and Jachc3, and was exploited by Girvin, MacDon- 
ald and Platzman in the magneto-roton theory of collective excitations of the incompressible 
states responsible for the fractional quantum Hall effected. Here we briefly discuss the part 
of this formalism that is most relevant to the present paper. 

We first consider the one-body case and choose the symmetric gauge. The single-particle 
eigenfunctions of kinetic energy and angular momentum in the LLL arec3'E3 

0m(^) = 7 — T^r^"' exp -^ , (17) 

where m is a non-negative integer, and z = {x + iy)/i- From (0) it is clear that any wave 
function in the LLL can be written in the form 

^(z) = f{z) e-^ (18) 

where f{z) is an analytic function of z^o the subspace in the LLL is isomorphic to the 
Hilbert space of analytic functionsEj'Ej'Ej. Following BargmannBa'El, we define the inner 
product of two analytic functions as 

{f,g)=Jdfi{z)r{z)g{z), (19) 

where 

dn{z) = {2ny^ dxdy e~^. (20) 

Now we can define bosonic ladder operators that connect ipm to (f)m±i (and which act on 
the polynomial part of (pm only) : 

at = -^, (21a) 

a = V2-^, (21b) 

oz 



so that 

a^ ipm = Vm + l(pm+i, (22a) 

a ipm = y/m ifm-i, (22b) 

if,a^g) = iaf,g), (22c) 

{f,ag) = {a^ f,g). (22d) 

All operators that have non-zero matrix elements only within the LLL can be expressed in 
terms of a and a^ . It is essential to notice that the adjoint of a^ is not z* j \/2 but a = \/2d/dz^ 
because z* connects states in the LLL to higher Landau levels. Actually a is the projection 
of z* / \/2 onto the LLL as seen clearly in the following expression: 

(/, 77f ^) = ("Tf ^' 3) = («^ /. 9) = (/, a g). 

So we find 

?=2|, (23) 

where the overbar indicates projecton onto the LLL. Since z* and z do not commute, when 
we need to project an operator which is a combination of z* and z, we must first normal 
order z*'s to the left of 2;'s, and then replace z* by z* . With this rule in mind and (p^), we 
can easily project onto the LLL any operator that involves space coordinates only. 
For example, the one-body density operator in momentum space is 

n — p-«q-r — p- 5 (9' 2+92*) _ -*- p-|g^* p-|gz 



where A is the area of the system, and q = Qx + i<ly Hence 



o~= — — e"*"^^ e"^^*'' = — — e~ 4 r ('24') 



where 

8 



r. 



9 



e-*'?ei-29 z ^25) 



is a unitary operator satisfying the closed Lie algebra 



r, 



rfc = r,+fce5^^^ (26a) 



[Tg, Tfc] = 2z Tg+fc sin -^r—, (26b) 



where g A A; = £^(q x k) ■ z. We also have TqTk T^qT^k = e^'^'^^ . This is a famihar feature of 
the group of translations in a magnetic field, because g A fc is exactly the phase generated by 
the flux in the parallelogram generated by q£^ and k£^. Hence the r's form a representation 
of the magnetic translation group [see Fig.(^]. In fact Tq translates the particle a distance 
£^z X q. This means that different wave vector components of the charge density do not 
commute. It is from here that non-trivial dynamics arises even though the kinetic energy is 
totally quenched in the LLL subspace. 

This formalism is readily generalized to the case of many particles with spin, as we will 
show next. In a system with area A and N particles the projected charge and spin density 
operators are 

p^= -^ ^e-q-» = _ ^e-^ Tq{{) (27a) 



1 A^ 1 N 2 

^." = ^ E e-'^- ^f = 7^ E ^"^ ^.(0 5f , (27b) 

where Tq{i) is the magnetic translation operator for the ith. particle and S'f is the /^th 
component of the spin operator for the zth particle. We immediately find that unlike the 
unprojected operators, the projected spin and charge density operators do not commute: 

|,-,.5Jl = iLe^^^2^^^SS;sm(^)/0. (28) 



This implies that within the LLL, the dynamics of spin and charge are entangled, i.e., when 
you rotate spin, charge gets moved. As a consequence of that, spin textures carry charged, 
as we will soon see. 



B. Spin-charge relation within the lowest Landau level 

We have argued above that in the SU(2) invariant case (i.e., when d/d = 0), the ground 
state at i/ = 1/m has its pseudospin fully polarized spontaneously. At these filling factors 
the interaction energy is minimized in such a state. For the same reason, we might expect 
that the low-lying excited states will be spin textures in which the local spin alignment 
varies slowly with position. To be explicit, we define the following as a spin texture state: 

|^[m(r)]) = e-^^|^o). (29) 

Here \iPq) is the S^ = N/2 member of the ground state spin-multiplet given in Eq.(0) and 
the operator O is a non-uniform spin rotation which reorients the local spin direction from 
z to m(r) (m is a unit vector). We limit ourselves to small tilts away from the z direction 

so that 

N 2 

= Ef^(r,)-S,=Ee^^^^~.' (30) 



where fl{r) = z x m(r) is the angle over which a spin is rotated. [Note that Q^{r) = 0, 
fi^(r) = — m^(r), and fi^(r) = ?7i^(r)]. We will later argue that our final result requires 
only that ^2 is slowly varying in space and not that Q is small. However the assumption 
is convenient at the present since it allows us to use a simple expression for Q and also to 
expand O as a small quantity. Projecting O onto the LLL ensures that \tp) has no projection 
on higher Landau levels^ as required in the strong perpendicular magnetic field limit. The 

|2 ^^ 

extra factor e^ in Eq.(^Ot) implies a non-standard definition for the Fourier components of 
Q^{r) which is adopted as a convenience. 

We can now calculate the excess charge density in a spin texture state: 

Spk = {He"^ Pk e-^^l^o) - {Hpkl^o). (31) 

Expanding in powers of O gives 

5pk = I {M[O.Pk]\i^o) - \ mP. [0,p,]]|^o) + . . . . (32) 

It is easy to check that the first term is zero. Using (p8|) we obtain 

[C, Pfc] = -7^ 2^ e 4 fi^ S^_q sm ^— , (33) 



substituting this into (|3^), and keeping only the second-order term we obtain 

5pu = -^ Ee ""'""'^""^ ^',K ^^^^{H[S%,WMo). (34) 

At this point we use our assumption that I'i/'o) is a state with uniform spin density so that 
the expectation value in Eq.(P^ is non-zero only when p = k — q. Using 

2 Jv 

[5'q_;j,5'^_J = — — ■ e 2 efj_ux2_^Sj, (35) 

i=i 

we obtain (using the fact that only the z component of spin is non-zero) 

^Pf^ = -^Jv^ ^M^ e-— E ^^ ^U sm -^. (36) 

Since fi is a slowly varying function, fi^ is negligible when q is large, we can make the 
approximation 

_ifcli . kAq kAq ,^_. 

e 4 sm ^~ -. (37) 

2 2 ^ ^ 

Substituting into (|36| ) we obtain 

= -^ ^ e,. Y.i^q fi^) A (^(A: - q) fi^J 
1 u 



'^ ^ ^M. E(vfi^), A {vn^k- 



~e^A{vn^)A{vnn\ . (38) 

10 



Here we have used the fact that N/A = z//(27r). Fourier transforming back to real space, we 
obtain 

6pir) = -^ e,, (vQ^ir) A Vfi'^(r)). (39) 

Expressing 5p in terms of m instead of ft, we finally obtain 

5p{v) = --—e m(r) • [(9^,m(r) x (9,,m(r)] (40) 

OTT 

which is exactly u times the Pontryagian index density, or topological charge density.t3'Eil In 



Eq.(40) we have used the fact that fi is small to replace z by m. The final result depends 
on VJl instead of Ct and it is clearS that the expansion in ( P^D is actually an expansion in 
terms of Vil, rather than ^2. Hence our final result is valid as long as ft is slowly varying 
so that Vri is small (compared to £~^).I13'EI 

The density in Eq. ( P0| ) can be viewed as the time-like component of a conserved (diver- 
genceless) topological 'three-current' 

ja ^ _^ ^a/37 ,^^^ m'^{r)dpm\r)d^m%r). (41) 

OTT 

Using the fact that m is a unit vector, it is straightforward to verify that d^^j^ = 0. 

Thus we have shown that for spin-states with S = N/2 the physical charge density is 
u times the the topological charge density, in the long wavelength limit. This remarkable 
result was first obtained by Sondhi et oL within the context of a Chern-Simons effective 
field theory description of spin textures.Ej The present derivation gives a microscopic proof 
of their result. The total extra charge carried by the spin texture is exactly the Pontryagian 
index: 

AN = -^ fdh t^, m(r) ■ [9^m(r) x 9,m(r)] , (42) 

which is an integer multiple of v because it is the number of times a unit sphere is wrapped 
around by the order parameter, i.e., it is the winding number of the spin texture.u For 
V = 1/m and m = 3, 5 elementary snin-textures carry the same fractional charge as the 
quasiparticles discovered by LaughlinEJ for spinless electrons. As we discuss below the fact 
that the charges are the same follows from very general considerations. Actually the spin 
texture states we have defined must contain precisely the same number of particles as {ipo) 
since the spin-rotation operator does not change the total electron number. However the 
spin-density may contain a number of well-separated textures with well-defined non-zero 
topological charge densities and hence well localized charges; only the net charge in the 
spin-texture states defined above will be zero. The system clearly also has states with 
locally non-zero net charge in the spin textures. 

The fact that spin textures carry charge can also be understood from the following 
very different point of view. In the Hartree-Fock picture the electrons see a very strong 
exchange field which locally aligns the spins ferromagnetically to produce the spin texture. 
As an electron propagates through this (slowly spatially varying) exchange field, its spin 
adiabatically follows the local orientation of the exchange field. A consequence of the tilting 



11 



of the spin is that when an electron moves along a closed path that surrounds the area F, 
the spin contributes a Berry's phase to the path integral^: 



V 



-- / d'^rtf,^ m(r) • [9^m(r) x 9^m(r)] . (43) 



This extra phase is exactly equivalent to having an Aharonov-Bohm phase due to additional 
magnetic flux inside F. In our system, the Hall conductance is not zero: 

^xy = -^, (44) 

which means that additional flux $ gives additional charge Q = eu^/^Q.M This is the same 
mechanism that causes Laughlin quasiparticles to carry quantized fractional charge when 
axy is quantized.^ Combining (^31) and (^4]) tells us that the additional charge density is 
given by Eq.(|^).Ea 



IV. EFFECTIVE ACTION 

A. SU(2)-invariant interactions 

In this section we calculate the effective action of a smooth spin texture for SU{2)- 
invariant (i.e. d/i = 0) electron-electron interactions. The considerations in this section 
apply to both a single-layer with zero Zeeman energy and to a pseudospin-polarized double- 
layer system in which the two layers are spatially coincident so that the interactions between 
layers and within layers are the same. In the double-layer case this limit can of course never 
be achieved experimentally, but it is a convenient place to begin the analysis since the ground 
state is exactly soluble in this limit. We assume, for the sake of convenience, that the spins 
are almost aligned in the z direction, and that they vary slowly in space, i.e., Qg is negligible 
when qi > 1. The interaction, after projection onto the LLL, is 

V = l EK(P.P-.-A^e-^), (45) 

where Vk = J d'^rV{r)e~^^'^. The expectation value of the energy is 

6E={^o\e^F,e-'^]\^Po) 

= -i{MV, 0]|^o) - I (MO, [O, F]]|Vo) + ■ ■ ■ (46) 

Since Q^ = 0, the leading term vanishes. The second-order term gives 

= -]Y.Vk {H\0, [O, p,p„J]|V^o). (47) 

^ k 

The commutators in Eq . (^7]) can be evaluated using 
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1 IfcP ^-^r , .. xr^ „,. — M 



[p„o] = ^e-^ j:[nij),j:n^^S' 



q --qi 



Q 



= ^iT.e''^^',S^r,_,{j)sin'-^. (48) 

Substituting this into the above expression, we obtain for small q 

k <? 

= I T^K^qK i-^<iW-q + {^q)^l {-^qW~q] 

^ q 
= Pl I £r [(Vfi")2 + iy^yf] = ^ I (friyuif . (49) 

The spin stiffness p°, implicitly defined above, is related to the pair correlation function of 

iV^o) by 

P's = ^2 j dkk'VkKk) (50) 

where the pair-correlation function h{k) = (z//27r) J d^r{g{r) — 1) exp(— ik • r). For u = 1, 
lipo) is known analytically and the pair correlation function can be evaluated analytically; 
h{k) = -exp(-|fc|V2). 

In this calculation we have kept the lowest order gradient terms only. The physical origin 
of the stiffness is the loss of exchange and correlation energy when the spin orientation varies 
with position. For the Coulomb interaction, p° = e"^ / {lQ\/27re£) ~ (e^/e£)2.49 x 10~^ at 
z/ = 1. For z/ = 1/3 and z/ = 1/5 we have evaluated p° numerically using hypernetted- 
chain-approximationcJ correlation functions and find p^ = (e^/ei) 9.23 x 10~^ and ps = 
(e^/ei) 2.34 x 10~^ respectively. For u = 1 this is exactlv the coefficient of the gradient 
term which Sondhi et alM obtained by fitting the knownE^'y long wavelength spin wave 
spectrum, but here we obtain it from a first-principles calculation. The classical model 
defined by Eq.(|i9|) is called the 0(3) non-linear sigma model and has been studied in great 
detail.Ea We note in passing that for the SU(2) invariant case, the spin stiffness ps found 
here is exact. Quantum fluctuation corrections to the Hartree-Fock approximation affect 
only higher gradient terms in the action. 

Now let us see what happens at higher order in Q. The third term is 

^^'^ = Y^y^h^MP, [o, [o,p,p_,]]]i^o). (51) 

k 

We have an odd number of powers of S^ or 5*^ combined together, but the state is polarized 
in the z direction, so the expectation value must be zero. In general odd order terms are 
zero by symmetry. The next non-zero term appears at fourth order: 

E^"^ = ^E^4^^°l[^' [^' [^' ^^P'^P-M^o). (52) 

k 

The nested commutators can be expanded as 
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[O, [O, [O, [O.VkP-M = PkPjOjO, [0,p_,,]]]] + 4[0, p,][0, [O, [0,J>_,]]] 

+Q\0,\0,p,]]\0, [0,p_fc]] +4[0, [O, [0,p,]]][0,p_,] 

+ [0,[0,[0,[0,pMP-k- (53) 

We remark that for z/ = 1 the first and last terms (with k y^ 0) give no contribution 
to E^'^\ since the ground state is annihilated by pk- (There are no nonuniform density, 
5*2 = N/2 states since we have a full Landau level.) The third order term contains the 
Hartree-interaction between the charges of the spin-textures as can be recognized after the 
following manipulations: 

WmMO, [0,p,]][0, [0,p_J]|Vo) ^IEMMO, [0,p,]]|^o)(^o|[0, [0,p_,]]|^o) 

^ k ^ k 

= 7jY.^k{Spk){Sp-k) 

^ k 



]-( fd^rd\'V{r-r')6p{r)6p{r'), (54) 



where dp = -^eij^diin x djin) ■ m is the topological charge density. In the first step above 
we have used the fact that IV'o) is an approximate eigenstate of [O, [0,p_k]], since fi^ is zero 
and even combinations of the spin operators S^ and 5*^ (i.e. S^ S^, S^ S^, etc) commute 
with S^. (For u = 1 this maneuver is exact.) Thus the third term in Eq. (|53D contains the 
direct Coulomb energy term. If we had a short range interaction, other terms would give 
contributions of the same importance, however because of the long-range of the Coulomb 
interaction, i.e., because V^ ~ 1/k as k ^ 0, one can show that this term is the dominant 
fourth order term in the spin-texture energy. We take this to be the next leading term in 
the spin-texture energy functional. 

So far we have calculated the most relevant terms of the static energy of a spin texture. 
The dynamics can be obtained by studying the equation of motion. The quantum equation 
of motion is 

^^ ^ {MO,[0,V]]m = ^,Jj^E[m], (55) 



where E[m\ is the energy functional of the spin texture, 

E[m] = ^ fd^riVmf + - f f d^d^r' V{r - r')6p{r)6p{r'). (56) 

If we include only the leading gradient term in the energy functional, an approximation 
which is always valid for sufficiently slowly varying spin-textures, we obtain 

o?mq 47rp^„ 

— , = — 1^ z X mq (57) 

dt hu ^ ^ ^ 

The equation of motion has spin- wave solutions in which the magnetization precesses around 
the z direction with wavevector q and frequency huj = {4tt p^q"^ / u) . This is precisely the 
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energy of the long-wavelength spin-waves of the system.tj'Ej. This equation of motion im- 
mediately leads to the following effective Lagrangian: 

L = — f rfVA[m(r)] ■ dtm{r) - E[m], (58) 

where A is the vector potential of a unit magnetic monopoleE§El in the spin space; i.e., 
Vm X A = m. [For spins oriented close to the z direction as in the above discussions 
A ^ (1/2) {—my, 771^,0)]. The first term simply contributes to the action a geometric phase 
proportional to the solid angle traced out by the spin vector during its motion. This is 
exactly the Berry's phase for the spin and appears at the adiabatic level as expectedE3. 

We have now derived all the terms in the effective Lagrangian of Sondhi et alxj from a 
first-principles calculation. 

B. Effective action for symmetry-breaking interaction 

In this section we derive an effective action suitable for double-layer systems using the 
pseudospin analogy discussed in Section |I|. We assumeCa that terms in the Hamiltonian 
where electrons are scattered from layer to layer by interactions can be neglected and that 
the two wells are identical. We define 

V'^l (V^.^ + V,^) (59) 



V: = I {Vjf - V,^) (60) 

where V/f^ is the Fourier transform with respect to the planar coordinate of the (intralayer) 
interaction potential between a pair of electrons in the same layer and Vj^ is Fourier trans- 
form of the (interlayer) interaction potential between a pair of electrons in opposite layers. 
If we neglect the finite thickness!^ of the layers, Vjf = 27ie'^/k and V/F = exp{—kd)Vif'. The 
interaction Hamiltonian can then be separated into a pseudospin-independent part with 
interaction V^ and a pseudosp in-dependent part. The pseudospin dependent term in the 
Hamiltonian is 

v,^ = 2j2v:sir,. (61) 

k 

Since V,f- > V/f, this term produces an easy-plane rather than an Ising anisotropy. The 
pseudospin symmetry of the Hamiltonian is reduced from SU{2) to ^(1) by this term. 
In addition, this term changes the quantum fluctuations in the system since it does not 
commute with the order parameter 

Fsb, 51 7^ 0, (62) 

where fi = x,y. 

The pseudospin texture energy due to the pseudospin-independent term in the Hamil- 
tonian can be calculated as discussed in the previous section. In this section we calculate 
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the contribution of the pseudospin-dependent term in the Hamiltonian to the spin-texture 
energy. 

In calculating the contribution of the pseudospin-dependent term to the spin-texture 
energy we approximate the ground state by the S = N/2 total pseudospin eigenstate which 
is the ground state in the limit of zero layer separation. We argue that the form of the 
energy-functional we derive must remain valid even when quantum fluctuations due to the 
pseudospin-dependent terms in the Hamiltonian are present. However, the coefficients which 
appear in the energy-functional will be altered by quantum fluctuations and the explicit 
expressions we derive below are accurate only when the pseudospin-dependent interactions 
are weak, i.e., only when the layers are close together. In Section |V111| we will discuss 
estimates obtained for the quantum fluctuation corrections to these coefficients from finite- 
size exact diagonalization and many-body perturbation theory calculations. 

It is convenient here to take the ground state {ipo) to be spin polarized along the x 
direction. We first calculate the energy change associated with small oscillations of the spin- 
texture away from the x direction. The leading term vanishes just as in the SU(2) invariant 
case. To understand qualitatively the physics contained in the pseudospin-dependent term 
in the Hamiltonian we focus first on the second-order term 

6E^^ = -J2V,^^o\2[0,Sl][0,S'.,] + S'.dO, PX]] + P, [0,r,]]Sl\^o). (63) 

k 

The first term on the right-hand side of Eq.(^) yields exactly (for u = 1) 

k k 



2^^ . 



^ \4^e 2 rnlm'-_ 



k"''-k 
k 

N^ 



2^2 



I d^rd^r'V' {r-r')m' (r)m" (r ') , (64) 



where V{r) is the Fourier transform of V^e 2". This is exactly the Hartree-like charging 
energy. In the limit of a smooth spin texture, using the gradient expansion, this term 
becomes a local mass term, which is the capacitive charging energy. That is. 



13h I d^r{m'Y (65) 

where 



2 
Ph = ^ I d'rVir). (66) 

We see immediately from this term that the symmetry-breaking interactions favor equal 
population of the two layers, or in pseudospin language they favor spin-textures where the 
pseudospin-orientation is in the x — y plane. 

The right-hand side of Eq.(^) can, after a straight-forward but lengthy and tedious 
calculation, be expressed in terms of the two-point correlations of iV'o)- The calculations are 
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similar to those in Ref . ( |^ and involve commutators of magnetic translation operators and 
of pseudospin operators. The following identity enters the calculation at several points: 

K^ Sg] = ^CA^. COS ^ e^P^ 5p+^ + ^-^ 6^, sm ^ e^^^ p^^^. (67) 

In the limit of slowly varying spin-textures we obtain the following result for the contribution 
of the symmetry-breaking term to the energy of the spin-texture 

E,b[m] ^ |rfV|/?„K)2 + ^(Vm^)2 - |[(Vm^f + (Vm^)2]|, (68) 



where 



and Pm = Ph + Pxc with 



7/ /"OO 



00 



Pxc = ^ I dkkV^h{k). (70) 

OTT"^ JO 

The total mass term is given by 

I3^ = (3h + I3^, = —J dk[V\Q)-V\k)]kh{k). (71) 

OTT Jo 

Notice that Va > Ve but that the intralayer interaction contains an exchange term which 
reduces the effect of Va- Also note that in the limit where the symmetry breaking interaction 
is local, i.e., the limit where Vj^ is independent of k, the exchange-correlation contribution to 
the coefficient of the mass term vanishes. We can understand this result in the pseudospin 
language by noting that in this limit the symmetry breaking interaction is proportional to 
J2i{SiY which (for spin-|) still commutes with all components of the total spin operator 
and does not, despite appearances, destroy the SU(2) symmetry of the Hamiltonian. From 
another point of view we can understand this result by noting that in lipo) no two particles 
can be at the same position and therefore they will not experience a local interaction. This 
property of iV'o) and the fact that jSm vanishes for local V;^ can be confirmed from the 
following identity satisfied by h{k): 

/■OO 

/ dkkh{k) = -V. in) 

h 

It is important to observe that in the limit of small layer separations V^f approaches vre^d 
which is local. Because there is no contribution to (5rn from this local term the mass coefficient 
ends up being proportional to the d^ rather than proportional to d at small d as would be 
expected based on naive considerations of the capacitance energy. 

Including both SU(2) invariant contribution defined in Eq.(|49|) and the symmetry break- 
ing contributions found above, the total energy-functional for a spin-texture is given by: 

ETotH ^ Jd\i^(3m{mn' + Y(Vm^)2 ^ Y[(Vm^)2 ^ (ym^)^]} (73) 
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where 



and 



PA = iT^ dkV,^h{k)k\ (74) 

6Iti Jo 



oo 



PE = ^2l dkV,''hik)k' (75) 

This result is easy to understand. The contribution to the exchange-correlation energy which 
is dependent on the z polarization of the pseudospin includes the Hartree energy which 
favors rriz = and the exchange-correlation within the layers. The exchange correlation 
energy within each layer increases superlinearly ~ p^/^ with the layer density so this term 
favors rriz 7^ 0. As discussed above the Hartree energy is always larger for constant spin- 
densities. The term proportional to (Vm^)^ in the energy density captures the reduction 
of the exchange-correlation energy from within each layer when the density in the layer is 
not constant and therefore pa is dependent only on the intra-layer interaction. [pA = p° at 
all layer separations. Because of the presence of the mass term, this gradient term is not 
important at long wavelengths.] On the other hand, pseudospin-order in the x — y plane 
represents interlayer phase coherence. As discussed earlier, an interlayer phase relationship 
which changes as a function of position results in a loss of interlayer correlation energy 
so that pe depends only on inter-layer interactions. In Fig. [§] and Fig. |^ we illustrate 
the dependence of P and pe on layer separation calculated from the above expressions for 
u = 1 and z/ = 1/3 respectively. We emphasize that these results are not expected to be 
accurate at large layer separations. We will compare these results with estimates from exact 



diagonalization calculations in Section VIII. 



C. Hartree- Fock Picture of Spin Textures and Gradient Expansion of Energy 

Functional 

In this section we develop a Hartree-Fock picture to describe spin textures and derive 
the corresponding energy functional. We also show that a gradient expansion of this energy 
functional gives exactly the result we obtained in the previous section. 

We work in the Landau gauge: A = (0, Bx, 0). The one-body orbital wave functions 
in the LLL in this gauge are: 

^x (r) = , ^ e'^yy e- ^^ , (76) 

where X = kyi"^ is the guiding center. 

A particular class of single Slater determinants at z/ = 1 in this gauge can be written in 
the form 

1^) = I[{ch cos ^ + Cl,, sin ^ e-(-)) |0), (77) 

where |0) is the fermion vacuum, C]^^ , creates an electron in the upper (lower) layer in orbit 
ipx respectively. In this state each Landau gauge orbital is occupied by a single electron 



whose pseudospin orientation is specified by the polar angles 6{X) and 0(X). Each Landau 
gauge orbital is localized within ~ £ of its guiding center. We are interested in states for 
which 9{X) and 0(X) vary slowly on the magnetic length energy scale so that \il)) describes 
a spin texture in which 

m^(x) = cos6'(a;), 

i7f{x) = sin^(x) cosip{x), 

m^{x) = sin^(x) sin(f{x). (78) 

lip) is not the most general spin texture, because m does not depend on y. As a consequence 
there is no spatial variation in charge density. 

It is straightforward to evaluate the energy oi\tp). For the following discussion we include 
the term in the Hamiltonian which allows electrons to tunnel from layer to layer and whose 
consequences will be explored in detail in a subsequent paper.EJ The Hamiltonian in this 
representation is 

H = f + V, 

T = -t^i CxiCxi + CxiCxi j , 
X ^ ^ 

^ = 2 2J ^X^X'iXiX^ ^XiCTi C'xaaa ^X^a2 C'xgai, (79) 

where a =t, i is the layer index, 

0T = yii = yA ^ y\i ^ yU ^ yE^ 

VX1X2X3X4 = j dridr2 1/(ri - r2)ip*x^iri)ipx^{r2)i^x3iT^i)i^xA^2)- (80) 

Hence we have 

(^|f| Vi) = -2t Y. COS ^^ sin ^^ cos if{X) (81) 

X 2 2 

= -t^m^(X) (82) 

X 



^y 



dx rrf{x) 



27r£2 

In the last step we allow the possibility that rrf depends on y as well. We see this result 
agrees exactly with the tunneling energy we obtained previously. The contribution from 
interactions can also be evaluated using Wick's theorem: 

(^IV^I^) = \ E VxlZx.x, Uc\^.,Cx,.,){C\^.,Cx,.,) - {C^x,.,Cx,.,){C\^.,Cx,.A 



^1-^2 ^3^4 



\ E I -^^(^1-^2) 
^ X1X2 1 
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+ [D^{Xi - X2) - Z^^(Xi - X2) - E^{Xi - X2)] m\Xi)m\X2) 

-E^(Xi - X2)K(Xi)m^'(^2) + m^(Xi)m^(X2)]| (84) 

In this equation we have absorbed the Hartree energy of the system for equal layer densities 
into the zero of energy. The quantities D{X) = Vx+yx,x+y,y and E{X) = Vx+y,yx,x+y are 
the direct and exchange two-body integrals for both intra and inter layer interactions. The 
above equation has clear physical content. The first term in the final form of the equation 
is the exchange energy in the absence of pseudospin polarization. The second term is the 
Hartree charging energy including a exchange correction. The third term is the exchange 
energy due to interlayer coherence, which is the source of the loss of exchange energy when 
m^ and m^ change. We now make a gradient expansion by writing 

m^{x 1)171" {X2) = [m'^{xi)]'^ + m'^(xi)(x2 — Xi) ( — m^(xi) j 

+ - m''{xi){x2 - xif {-K-^m''{xi)j H , (85) 

Replacing the summation over guiding centers by integrations, we easily obtain recover 
Eq. [ ^ . (The coefficients of the gradient terms are proportional to the second moments of 
the exchange integrals, J2x ^"^ElX).) In this formulation we see explicitly that the leading 
gradient corrections are adequate as long as the pseudospin orientation changes slowly on 
the scale of £. 

D. Hubbard-Stratonovich transformation approach to the effective action 

In subsection B above, we derived the effective action by calculating the energy functional 
for spin textures. The dynamical term in the effective action was determined by requiring 
that the Lagrangian implies the correct equation of motion for the spin textures. In this 
section we briefly sketch an alternate route for deriving the same effective action. The basic 
idea is the familiar Hubbard-Stratonovich (HS) transformation. We introduce auxiliary 
fields to decouple the interaction and integrate out the fermionic degrees of freedom to 
obtain the effective action of the auxiliary fields. The auxiliary fields are essentially the 
order parameters. 

It is inconvenient to project onto the lowest Landau level until the end of the calculation 
so we work with the full Hamiltonian which has the form 



H = Ho + V, (86) 



where 



^o = E7r- Pj--M^jJ 



J 2m 



e 
c 



(87) 



is the kinetic energy of the particles in the presence of the external magnetic field. This 
approach is more easily implemented if we have a delta- function-like interaction V{r) = 
Vo6{r) and so we discuss this case first: 
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y =^Y. PlP-1 = -^0 E Sg ■ S^g - ^ ^ PgP-q. (88) 

^ g q ^ q 

The advantage of a delta function interaction is that it can be expressed in terms of spin 
operators,E2l which makes it possible to decouple the interaction in terms of spin auxiliary 
fields. Written in the above particular form, we will find that the saddle point of the auxilian^ 
fields corresponds to the Hartree-Fock mean-field Hamiltonian.ell The partition function iglHl 

Z(/?) = Tr{e-^^°^(/?)}, (89) 

where 

5(/3)=T,e-/o''^M, (90) 

and 

V{t) = e^'^Ve-^°\ (91) 

Tr is the (imaginary) time ordering operator. Now we can introduce a vector auxiliary field 
h(r, r) and a scalar auxiliary field 0(r, r) to decouple V: 

After the HS decoupling, we find the fermionic Hamiltonian becomes quadratic so we 
can (at least in principle), carry out the trace over the fermion degrees of freedom and hence 
obtain the effective action in terms of the auxiliary fields h and 0. In doing that, however, 
we still need to make approximations. We notice that the direction fluctuations of h are 
massless (due to broken SU{2) symmetry) while the fluctuations of are massive.El It is 
therefore a good approximation to assume that (p is "frozen" to be a constant in space and 
time so it only contributes a chemical potential like term in the fermion Hamiltonian and is 
hence unimportant. Thus we will concentrate on the fluctuations of h. 

In computing the trace we notice that the trace is over nothing but the propagator of a 
system governed by the time dependent Hamiltonian 

H(t) =Ho- f (fr S(r) ■ h(r, r) + const. (93) 

in imaginary time. Since there is always a large one body gap for H{t) while the collective 
modes of h are gapless (which means H{t) is slowly varying), we can use the adiabatic (or 
Born-Oppenheimer) approximation to evaluate the propagatorE3: 

Tr{e-^^M} = e^^[^l-/o'^W'^^ (94) 

where 7[r] is the Berry's phase determined byE3 

z^ = -(vl/(r)||-vl/(r)), (95) 
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where |^I/(t)) is the ground state of H{t) and E{t) is the ground state energy. In order to 
proceed, we need to find |^(t)) and E{t). It is now possible to take the strong magnetic 
field limit. It is both plausible and easily checked that in this case |^1'(t)) is nothing but our 
familiar spin texture state: 



|v[/(r)) = e^oM|M/o), (96) 

where |\E'o) is the fully polarized state in z direction and 

0(r)=^S,.f^(r,rj), (97) 

j 

where 0(r, r) = z x m(r, r) and ni(r, r) is a unit vector in the direction of h. 
Now let us look at the Berry's phase term: 

id'y = 1- (^(r)l^(r + dr)) 

= 1 - (^(r)|e*^(^+'^")e-*^(^)|^)(r) 

= 1 - (^(r)|e^^^^(^)'^^"5[0(r),9.0(r)]dr|^(^^^)^ 

= -z{^{r)\dMr) + '-[0{T),drO{rm{r))dT 

= i^dT f d^r A{m{T,r)) ■ drm{T,r). (98) 

Here p = 1/{2ti(?) is the density of electrons. Hence we find that the Berry's phase term is 
exactly the dynamical (time dependent) term we obtained previously: 



^7 
Now let us calculate E(t): 



d f dr f dh A{in{T,r)) ■ drni{T,r). (99) 



E{r) = (M/(r)|iJ(r)|M/(r)) = (M/o|e-^«MiJ(r)e^«M|M/o) 
= {^o\H - t[0, H] - i[0, [O, H]] + --- \^o) 



fdhh-m+ ^pf I dh (V/i^) ■ {Vm") + ■■■. (100) 



The first term is just the gain of energy by polarizing the spins in the direction of the external 
field. The second term comes from the nonlocality induced by the projection to the LLL, so 
the spins cannot take full advantage of the external field if its direction is changing in space. 
This is exactly the physics that is responsible for the stiffness. Hence the effective action is 

S^l^]=i^[T]+ ["dr f dh-A' - ^h + lpf{Vh^) ■ (Vm^^)}. (101) 

Jo J 4 2 o 

Here h = |h|. We find the action has its minimum when 

h = pVoh, (102) 
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where n is a constant unit vector. At this point the one body Hamiltonian given by h is 
exactly the Hartree-Fock Hamiltonian. The fluctuations of the magnitude of h are massive 
and we can integrate them out to obtain the non-linear sigma model: 

SE[m]=t^[T] + ^PsfdTld^r{Vin)\ (103) 

where the stiffness ps = p^£^Vo/4 = Vo/(167r^£^), is exactly what we obtained from the spin 
texture calculations. 

The above discussion was for the case of short range interactions. To generalize to the 
finite-range-interaction case we separate a 6 function like part from the interaction, follow 
the same procedures as illustrated above, and treat the remaining part of the interaction 
as a perturbation. Note that non-delta function like interactions are not necessarily SU{2) 
invariant. At the first order in this perturbation theory (where {"^{t)) is unaffected), we 
recover all the results obtained in previous sections. 



V. CHARGED OBJECTS IN THE SYSTEM 
A. Skyrmions 

We start by briefly reviewing and commenting on results for charged excitations obtained 
by Sondhi et a/.t3 for a single-layer system of spin-1/2 particles with no Zeeman coupling. 
These results can be directly taken over to the double-layer system with d/i = (which is 
therefore SU{2) invariant). In this case one findscj'Ej, as discussed in the following section 
and mentioned previously, that the double-layer system with spontaneous interlayer phase 
coherence supports neutral gapless 'spin-wave' excitations which disperse quadratically in 
the long wavelength limit. This property of spin waves is characteristic of isotropic fer- 
romagnetism. However, in contrast to the case of quantum Heisenberg ferromagnets on a 
lattice, quantum Hall systems also possess charge degrees of freedom and are more analogous 
to itinerant electron ferromagnets. For example, after a spin is flipped, it can be scattered 
to other orbital states and carry charge throughout the system. In addition to the gapless 
spin wave excitations, there are gapful charged excitations. Only charged excitations can 
contribute to the conductivity so that the low temperature transport coefficients will be 
activated with an activation energy which is half the charge gap. Some of the low-energy 
charged excitations can be generated from the topologically non-trivial spin configurations 
discussed below. 

We have seen in Section HI that the physical charge density of the system aX v = 1/m 
is 1/m of the topological charge density. As a consequence of this the topological solitons, 
i.e., skyrmions, carry 1/m units of physical charge (see Fig. [^]).E3ll3 Inside a skyrmion the 
spins wrap around the unit order-parameter sphere exactly once.E3 In the SU{2) invariant 
case with Coulomb interactions, Sondhi et. al. have shown that skyrmions are the lowest 
energy charged excitations of the system at i/ = l.t3 This is not surprising because for the 
skyrmion spin configuration, the spins are nearly parallel locally, so the exchange energy is 
only slightly reduced. In contrast, for ordinary single particle excitations [see Fig.(P)], an 
added electron has its spin opposite to the others and has no exchange energy. As pointed 
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out earlier, in the SU(2) invariant limit we know the exact spin stiffness. Hence the exact 
energy of a single (large scale) skyrmion can be obtainedt3ll3: 

Es = 47rp,. (104) 

For the case of a system with Coulomb interactions at z/ = 1, we obtain from the non-linear 
sigma model energy expression 

It is important to realize that since the total particle number is fixed in our derivation 
of the non-linear sigma model energy expression, Eg actually gives the energy to introduce 
an isolated skyrmion or anti-skyrmion into the bulk of the system and maintain charge 
neutrality by adding or subtracting charge from the edge of the system. (For related careful 
discussions of quasiparticle energies in the fractional quantum Hall effect see Refs.[ ^,|^.) 
This energy must be subtracted off if we wish to calculate e±, the energy to add (+)or 
subtract (— ) 1/m electrons from the system aX v = 1/m in the form of a skyrmion or 
anti-skyrmion spin-texture. It follows^ that 

e± = ±z/e(z/) + Es (106) 

where ^{u) is the energy per electron in the incompressible ground state at z/ = 1/m. 
(For the Coulomb interaction case ^(1) = -^f7/8{e^/i) ^ -0.6266(eV£) audi ^(1/3) ^ 
— 0.4100(e^/£).) The chemical potential for u > 1/m is me+ while the chemical potential 
for u < 1/m. is —me-. For z/ = 1 it follows from the above that 6+ = — (l/4)i/7r/2(e^/£) 



and e_ = (3/4)W7r/2(e^/£). The energies of the localized quasiparticle excitations of Fig.0 
cannot be reliably calculated from the non-linear sigma model energy expression but for 
V = 1 the microscopic energy can be calculated analytically. For these excitations e^ = 
(the spin-reversed added electron has no exchange energy!) and e_ = J-n /2{e^ /€). In 
the absence of Zeeman coupling it follows that the lowest energy particles and holes are 
both formed from skyrmion spin textures. The lowest energy particle-hole excitation is a 
skyrmion anti-skyrmion pair which has energy lEg. This is only one half of that of the 
ordinary particle-hole pair in the case of the Coulomb interaction.E3 These results receive 
unequivocal support from numerical calculations.E3 Using our result for the spin-stiffness we 
can extend this analysis to the case of z/ = 1/3. At this filling factor Es ~ 0.0116(e^/£), 
e+ ~ — 0.1251(e^/£), e_ ~ 0.1483(e^/£). In the limit of large Zeeman energy at this filling 
factor both the localized quasihole and quasiparticle excitations will be completelvspin- 
polarized. The quasiparticle and quasihole energies in this limit have been estimatedaca to 
have the values e+ ~ — 0.120(e^/£) and e_ = 0.2337. Again the skyrmion and antiskyrmion 
states possible at zero Zeeman coupling have lower energy, although only barely so in the 
quasiparticle case. The particle-hole creation energy lE^ for v = 1/3 is approximately four 
times larger than in the large Zeeman coupling limit. This is consistent with results from 
the finite size exact-diagonalization study of Chakraborty et. al. who foundc3 that in the 
absence of Zeeman coupling, quasiparticle energies at z/ = 1/3 could be reduced by fiipping 
a single spin. 
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Another consequence of skyrmions being the lowest energy charged excitation is that in 
finite size systems on a sphere, the total spin of the ground state changes suddenly from y 
to zero or one half (for odd and even N respectively) when one changes the particle number 
from A^ to A^±l, where N is the Landau level degeneracy without spin.E3'E3 (We have verified 
that the same property holds for u = 1/3 in agreement with the analysis of the preceding 
paragraph.) This is because when a skyrmion is put on a sphere, the spin configuration is 
like a hedgehog and it is plausible that the total spin is essentially zero.Ej [As we discuss 
below the total angular momentum J is precisely zero.] Note however that the system is 
still ferromagnetic in the sense that it is not a local singletE3 and its Zeeman susceptibility 
still diverges. 

In this case the topological charge density is uniformly distributed on the sphere and 
the skyrmion is unfrustrated. The situation is very different for a system with rectangular 
geometry with periodic boundary conditions, i.e, the geometry of a torus. In this case the 
skyrmion is frustrated and hence squeezed by the finite size effect, so its size is much smaller 
than the system size. This can be understood by looking at the effect of periodic boundary 
conditions on the energy of a skyrmion. We know the gradient energy term is scale invariant, 
and for the ideal skyrmion solution, the energy is minimized to be 47rps- However an ideal 
solution does not satisfy the boundary conditions. In a rectangle a skyrmion has to be 
distorted near the boundary so the energy from the stiffness term will increase. It is obvious 
that the smaller the size of the skyrmion, the smaller the energy cost due to the boundary 
effect will be. However the skyrmion cannot be too small because it costs too much Coulomb 
energy (which wants the skyrmion to be as large as possible so that the excess charge will 
be distributed as uniformly as possible). The optimal size of a skyrmion on a rectangle is 
determined by a competition between stiffness and Coulomb energies and as a consequence 
the total spin S of the ground state of the A^ + 1 particle system is size dependent. 

We can estimate the difference between (A^ + l)/2 and S (i.e., the number of spins flipped 
in the ground state) in the following way. Keeping the two leading terms, the energy of a 
spin texture at z/ = 1 is 

E = ^p, f dh ( Vm)2 + 1 / (frdh' V{r - r')p(r)p(r') , (107) 

where p is the topological charge density and V is the Coulomb interaction. The energy of 
a skyrmion with linear size A 3> £ in a system with linear size R is 

Ei\)=A7rp, + A{^r + B^, (108) 

where A and B are positive constants with units of energy. The flrst term in Eq.( p.08|) is 
the usual energy of an unfrustrated, inflnite size skyrmion, and the last two are flnite size 



corrections from the two terms in Eq. [107]. respectively. Minimizing E with respect to A 
gives 



2 „1 



\(xm&. (109) 

The number of spins flipped in a skyrmion with size A and long distance cutoff R isi^l 

/\So^\^\n{^)o^Ni\nN, (110) 

A 
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where we used the fact that A^ oc R^. From Eq. ||110|| we see that as N becomes large, 



the number of flipped spins in the optimal skyrmion state gets large but is always small 
compared to A^ so that the ground state is almost fully polarized. 

Skyrmion spin-textures produce an excitation energy which is independent of the texture 
size as long as the size is large compare to microscopic lengths {£) and small compared to the 
system size. This leads to dramatic finite-size effects which are typified by the qualitative 
differences between ground state spin quantum numbers for electrons in a rectangle and on 
a sphere, where the state is unpolarized. For the case of a square with periodic boundary 

2 

conditions, a crude estimate gives the coefficient in front of N^lnN in Eq.(|ll(JD to be 



0.26. In Fig.(|^) we compare this estimate with results from finite-size exact diagonalization 
calculations for square boundary conditiongl^ and find qualitative agreement. Using the 
composite fermion theory, Jain and WuEl gave an alternative explanation of the fact that 
the total spin goes to zero for an A^ + 1 electron system on a sphere. Their theory does not 
distinguish between the sphere and torus and we do not believe the results on a torus can 
be easily understood in their formalism. The quantitative agreement with the ground state 
energy on the sphere and the qualitative agreement for the polarization on the torus lend 
strong weight to the skyrmion picture. 

We close this discussion by noting that it is possible to write down simple microscopic 
variational wave functions for the skyrmion, both in the plane and on the sphere. Consider 
the following state in the plane 



^A = n r ^v^' (111) 




where \l/y is defined in Eq.(|^), ()„ refers to the spinor for the mth particle, and A is a 
fixed length scale. This is a skyrmion because it has its spin purely down at the origin 
(where Zm = 0) and has spin purely ud at infinity (where Zm ^ A). The parameter A 
is simply the size scale of the skyrmiont3'E3. Notice that in the limit A — ^ (where the 
continuum effective action is invalid, but this microscopic wave function is still sensible) we 
recover a fully spin polarized filled Landau level with a charge-1 Laughlin quasihole at the 
origin. Hence the number of flipped spins interpolates continuously from zero to infinity as 
A increases. 

In order to analyze the skyrmion wave function in Eq.( |lll|) , we use the Laughlin plasma 
analogy. In this analogy the norm of tp\, Trjo-} / D [z] | \I' [z] p is viewed as the partition 
function of a Coulomb gas. In order to compute the density distribution we simply need to 
take a trace over the spin 

Z = /'/)[^]eK^»>^-'°^l^'-^^l+^^^'°^^l^^l'+''^~^^^l''=l'}. (112) 

This partition function describes the usual logarithmically interacting Coulomb gas with 
uniform background charge plus a spatially varying impurity back ground charge Aph{r), 

Ap.(r) . -lv^V(r) ^ -J^^J^. (113) 

V{r) = Uog{r^ + \^). (114) 
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For large enough scale size A ^ £, local neutrality of the plasmacll implies that the 
excess electron number density is precisely ApJr), so that Eq.( [114| ) is in agreement with 
the standard result for the topological density.El 

For a complete microscopic analytic solution valid for arbitrary A, we use the fact that 
the proposed manybody wave function is nothing but a Slater determinant of the single 
particle states 4>m{z), 

^27r2™+im!(m + 1 + ^) W 
The electron number density is then 

Af-l 1 

AnW(^)^5:i0„(z)p--, (116) 

m=0 ^^ 

which yields 

ZTT V2 JO / 

Similarly, the spin density distribution S''^{r^ can be obtained, and it also agrees with that 
for the standard skyrmion inE3 the large A limit. We have also computed the skyrmion 
creation energy from the spin-dependent pair correlation functions of the plasma following 
the same procedure as in Ref.( |^). Fig.(|^) shows a plot of this energy as a function of 
scale size A and shows that the microscopic formula gives the correct asymptotic value of 
one-half the quasi- hole energy for the large A-limit, in which the continuum field theoretic 
picture holds exactly. As the core size decreases, the skyrmion energy increases due to the 
increasing Coulomb charging energy. However it does not diverge as the naive extrapolation 
of the field theoretic expression would. 

Finally, we note that by replacing \\\ by (^„), we can generate a skyrmion with a 
Pontryagin index n. 

The skyrmion wave function has a particularly simple form on a sphere. On a sphere 
with radius R = S^^'^i where S is an integer or a half integer, the number of a flux quanta 
is Ng = 2S. The single particle kinetic energy ist3 

T = ^uj,\A\yS, (118) 

where 

A = r X (-iV + eA(r)) (119) 

is the kinetic angular momentum. One can show that L = A + Sfl (where Jl = -ct) is the 
generator of rotations for the system, i.e.: 

[L'',X^] = ieafs-yX\ (120) 

where X is any vector. We also have 
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\A\^=\L\^-S\ (121) 

so the eigenvalues of |Ap have the form {n + S){n + S +1) — S"^, where n is an integer. For 
n = 0, one obtains the LLL energy T = ^Uc, and the degeneracy is 25 + 1 = A^^ + 1. If we 
use the Dirac gauge A = ^0cot6', everything can be easily expressed in terms of spinor 
coordinates: u = cos | exp(i|), v = sin | exp(— z|). In this represent at ioncj 

L+ = ud^, (122) 

L- = vdu, (123) 

L' = ^iudu-vd,). (124) 

The LLL wave functions are simply homogeneous polynomials of u and v of degree 25*. The 
filled LLL single Slater determinant is just 

N 

Y[{uiVj-UjVi), (125) 

where N = Ng + l. The single antiskyrmion (that carries charge -1) wave function is simply 

Af-l / \ N-1 

fe=l V "-k J i<j 

where we have explicitly put in the fact that the total number of particles is now A^ — 1. 
The spin configuration of ipas is that of a hedgehog (with spins pointing inside toward the 
center of the sphere) since the ratio of \uk\ to \vk\ varies as cot(6'/2). This state is neither 
an eigenstate of Stot, nor an eigenstate of Ltot- It is however a singlet of the total angular 
momentum J: 

J = Ltot + Stof (127) 

Physically this means the state \ipas) is invariant under a spin rotation followed by an identi- 
cal space rotation, which is clear from the uniform nature of the hedgehog spin configuration. 

If we project xjjas onto the subspace of Stot = (or | if we have an odd number of 
particles), we automatically get Ltot = (or i). So the projected state will be invariant 
under both spin and space rotation. This state should have good overlap with the exact 
ground state. 

The skyrmion (that carries charge +1) wave function has a similar form: 

N+l /a \ N+1 

fc=i \ ^^k ) i^j 

The spin configuration of this state is exactly the opposite of '^as-, i-e., it is like a hedgehog 
with all spins pointing outward. 
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B. Merons 

When d/ a, 7^ 0, the z component of the order parameter is massive and the system 
has U{1) symmetpj. In this case, there is another class of topologically stable charged 
objects, meronsSMta As shown in the following, merons [see Fig.(P)] carry one half unit of 
topological charge and henceE^ l/2m units of electron charge. Far away from the core of a 
meron the order parameter lies in the (massless) XY plane and forms a vortex configuration 
with ± vorticity, while inside the core region the order parameter smoothly rotates either 
up or down out of the XY plane. Hence there are four flavors of merons. The energy of 
a single meron diverges logarithmically with the system size with a coefficient proportional 
to the inter-layer spin stiffness. The interaction between merons has a contribution from 
the stiffness energy which is also logarithmic, attractive for opposite vorticity pairs and 
repulsive for same vorticity pairs. These properties are exactly the same as the vortices in 
the classical XY model. In order to determine the sign of the charge carried by a meron, one 
has to specify both its vorticity and the spin configuration in the core region. Merons will 
also have a long range 1/r interaction due to their charges which is attractive for oppositely 
charged merons and repulsive for like-charged merons. 

The fact that merons carry topological charge one half can be seen by the following 
argument. Imagine a vortex in the spin system. If an electron circles the vortex at a large 
distance, its spin rotates through 27r. This induces a Berry's phase of exp(i27rS') = —1 which 
is equivalent to that induced by a charge moving around one-half of a fiux quantum. Since 
<^xy = e^ /fnh, the vortex picks up charge l/2m. The topological charge of a meron can also 
be understood by considering a variational function for the meron spin texture: 



m 



y 1 — {m'''{r)Y cos 99 , yl — (m^(r))2siny9 , m^{r) > . (129) 



The local topological charge density calculated from 6p = — ^ejj(5jm x djicn) ■ m can be 
expressed in the form 

5p{r) = -L^, (130) 

Anr ar 

and the total charge is 

Q= fd\6p{r) = -K(cx))-m"(0)]. (131) 

For a meron, the spin points up or down at the core center and tilts away from the z 
direction as the distance from the core center increases. Asymptotically it points purely 
radially in the x — y plane. Thus the topological charge is ±1 depending on the polarity of 
core spin. The variational function mentioned above corresponds to a vortex with positive 
vorticity. In order to make a vortex with negative vorticity (anti- vortex), we need to apply 
the space-inversion operation to the vortex solution. Since topological charge is a pseudo- 
scalar quantity, it is odd with respect to parity. Hence the general result for the topological 
charge of the four meron fiavors may be summarized by the following formula: 

Q = ^K(oo)-m^(0)]nv, (132) 
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where n^ is the vortex winding number. 

Finite energy excitations can be formed by pairs of merons with opposite vorticity. It 
seems hkely that under appropriate circumstances the lowest energy charged excitations of 
the system will consist of a bound pair of merons. (A skyrmion can be viewed as a closely 
bound pair of merons with the same charge and opposite vorticity and a meron can be 
viewed as half a skyrmion.) The energy of a pair of merons with opposite vorticity (but like 
charge) separated by a distance R is given by 



Emp = 2E^, + -rB + '^T^PE HR/Rmc) (133) 

where E^^ is the core energy of an isolated meron R^^ is the core radius of an isolated meron, 
and the expression should be applicable only when R ^ R^c- Minimizing this expression 
with respect to R gives a meron separation R* = e'^/SirpE- Using the expression for the pe 
derived in the last section this gives R* ~ 6i for d/i ~ 1. Quantum corrections are expected 
to reduce pe so this expression should give a lower bound on the meron separation. This very 
attractive picture of the lowest energy charge carriers in the system [see Fig.(|T^)] will only 
apply when the meron separation is larger than the meron core size (which is expected to be 
~ £), and its energy is lower than the energy of a conventional quasiparticle excitation where 
a charge is added with pseudospin directed in opposition to the local pseudospin order. It is 
clear that we should expect the meron core size to increase as d/i approaches zero and the 
mass term in the energy expression becomes small since there is only a small energy cost for 
pseudospins to point out of the xy plane. Hence the picture is not likely to apply for very 
small d/i. Further work which estimates meron core energies and radii will be necessary to 
substantiate this picture and is currently in progress. We note in passing that the above 
description of charged vortex antivortex pairs can also be used mutatis mutandis for neutral 
vortex antivortex pairs. These will have a conserved momentum and the neutral collective 
mode (discussed in Section |V|) will cross over from spin waves to such 'quasiexcitons' at 
large wave vectors in analogy to what occurs in the single layer case.c2l 

As in the case of skyrmions, we can write down explicit microscopic variational wave 
functions for vortices (merons). We start with the simplest example: a meron with vorticity 
+1 and charge — | that has the smallest possible core size: 

M ^ 1 

1^+1,-i) = n (yl^T + y|4+u)|0). (134) 

Here |0) is the fermion vacuum, cl^i- < creates an electron in the upper (lower) layer in the 
angular momentum m state in the LLL and M is the angular momentum quantum number 
corresponding to the edge. The vorticity is +1 because far away the spin wave function is 
essentially 

X(0) = ^ ( '7 ) , (135) 

where (p is the polar angle. The charge is — | because we have created a hole in the center 
of the lower layer (m = | is unoccupied). Since the spin is pointing up at the center, 
this agrees with the spin-charge relation derived earlier. From the spin-charge relation we 
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know we can change the sign of the charge of a meron by changing the direction of spins in 
the core region without changing the vorticity. This can be seen exphcitly from the wave 
function: 

M ^ 1 

|^+i,+ i) = cl n (^4t + ^4+u)|0). (136) 

This state has charge +| because we have put in an electron in the m = state in the lower 
layer. Obviously what we did is to flip the spins in the core region to the down direction 
without changing the vorticity of the meron, in the spin texture language. 

A general wave function that describes a meron with vorticity k has the form 

M 

|^+fc,_|) = n («™4t + bmcl+,^)\0). (137) 

This meron has charge — | and we have assumed k > 0. Generalization to other cases is 
trivial, a^ and b^ are parameters that satisfy 

I |2 I 11, |2 -I 

I ^m I ~r 1 0,jj I i 

lim ^ = e^-^o (138) 

m— >oo h 

where 0o is a constant. By adjusting the asymptotic behavior of a^ and bm (in particular 
their ratio, Eq. (|137|) can also describe skyrmions and other charged objects, while adjusting 



these coefficients in the core region one can modify the spin configuration there. Since ( |137| ) 
is a single Slater determinant, we can calculate its energy using the original Hamiltonian 
directly without using the effective energy functional derived earlier. After subtracting 
the energy of the fully spin-polarized state, we can show that the vortex energy diverges 
logarithmically with the system size as would be expected for an XY system: 

AEk^jk^lnM (139) 

where 

oo 9/ -I- 1 

7 = E i-^y-^vf- (140) 

1=0 ^ 

Here the vf are Haldane's pseudopotential parameters for the interlayer interactionE3, and 
we have neglected terms that are finite. If we use the energy functional, we find the divergent 
part of the vortex energy is 

AEk^^PEkHnM (141) 



implying 



oo 9/ -I- 1 

PE = E(-l)'^-f- (142) 



1=0 



One can easily show that this expression of pe agrees with previous expression exactly. The 
present expression in terms of pseudopotential parameters may be more useful in finite size 
numerical calculations.E3 
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VI. COLLECTIVE MODES AND PSEUDOSPIN RESPONSE FUNCTIONS 



In this section we combine the equation of motion for the spin-textures (Eq.|^3[) and 
the spin-texture energy functional (Eq. |7^) to calculate the pseudospin linear responseSH 



functions for d/d. ^ 0. We take the pseudospin of the system to be polarized in the x 
direction and consider the linear response to a time- and space-dependent Zeeman field in 
the y — z plane. Fourier transforming with respect to both time and space we find that 



-luo 



-^(2/5 + gVA)\ (my\_(-K 



where hy and h^ are the Fourier coefficients of the pseudospin magnetic field at frequency uj 
and wavevector q. Physically hz corresponds to a time- and space-dependent bias potential 
between the two wells, while hy could arise from a space- and time-dependent interlayer 
tunneling amplitude. We see immediately that the response is singular when 

/47r\ ^ r ^ 

u^ = ul^ = [—^ [213 + q^pA\ q^PE. (144) 

Here uocm is the frequency of a long- wavelength collective mode of the system. For the d/i, = 
case , P = 0, Pa = Pe = p'i, and the collective mode frequency reduces to the result obtained 
previously for the spin-wave collective mode of isotropic ferromagnets. [ujcm = 47rg^p°/z/.] 
The collective mode corresponds to a spin-precession whose ellipticity increases as the long- 
wavelength limit is approached. The presence of the mass term (/? ^ 0) changes the collective 
mode dispersion at long wavelengths from quadratic to linear. In the limit of small q 



V 



— J2(3pEq. (145) 



We can solve Eq. ||143|| for the frequency and wave vector dependent linear response to 



bias potential between the two wells. The result is 



A-nq^PElv ,^.ps 

Xzz{q,^) = ^ —■ (146) 



'cm 



It is interesting to compare this with formally exact relations for this response function 
which can be obtained from the microscopic Hamiltonian of double layer systems. We first 
note that Xz^ is related to the dynamic structure factor, 

SzMM = J2\{<iIn\S'g\'^o)\'S {hu - {Er, - Eo)) , (147) 

n 

by 

Szz{q, uj) = Imxzziq, uj + irj). (148) 

71 

Here |\E'„) is an exact eigenstate of the double-layer system. Our result for Xzz{q,^) thus 
implies that at long wavelengths 
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...(g, -) = 1^ (2^ f^.^^/ (^ - ^c.(g)) . (149) 

Some frequency moments of this structure factor can be related to ground state correlation 
functions of the double-layer system. The first moment gives the oscillator strength. From 
Eq. [p9f| we find that 

fzz{q,^)= dujujs^^{q,uj) = — q'^pE- (150) 

JO z/ 

This result agrees with results for this moment calculated previouslyt9'u directly from the 
microscopic Hamiltonian as we can confirm using Eq.||7^. We see that the oscillator strength 
vanishes like g^ as expected for an ordinary superfluid, however the coefficient is proportional 
to pe and hence is non-universal. [The remaining oscillator strength is found in a high 
frequency collective mode lying below omegac by an amount proportional to pe- as shown 
below in Eq.(^).] 

The zeroth moment of the dynamic structure factor gives the static structure factor: 



(g) ^ / dcosUq,^^) = {^o\S\S^J^o) = SWttt^T^^V^- (151) 




We note that the static structure factor vanishes linearly with q as q goes to zero. This 
property illustrates a qualitative differenceE3 between the ground state of the double-layer 
system at d/i ^ and the ground state in the dji = limit for which s^ziq,^) approaches 
a constant as q goes to zero. This property is analogous to what happens in a repulsively 
interacting Bose gas such as ^He in which the structure factor vanishes linearly with q. 
According to Feynman's single-mode approximation picture, this is required in order to 
achieve a linearly dispersing Goldstone mode.E2l The proper structure factor can be included 
to improve the variational ground state wave function by means of the Jastrow ansatz 



^ = exp - ^ -Sl^S^^ n L./2 ^v^, (152) 




where A is a variational parameter. The incompressible mmn states with n ^ m discussed 
in section fXl have the propertyEil that Szz{q) ~ q^ at long wavelengths. This property holds 



for the ground state whenever there is an excitation gap and additional Jastrow factors are 
not required to capture the correct long length scale fluctuations. 

For d/i, 7^ the mass term suppresses long length scale fluctuations in S^ which mea- 
sures the difference of the density in the two-layers. The minus-one moment of Szziq-,^) is 
proportional to the static response function: 

^-^^'" = °^ = M¥T7^ ^'''^ 

Note that Xzziq-,^ = 0) diverges in the d/i limit because of the broken SU{2) symmetry in 
the ground state of the d/i = system. 

The collective mode dispersion can also be obtained from a Lagrangian formulation which 
may be useful in describing the thermodynamics of the system. The partition function can 
be expressed by 
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Dme-^ [m] ^;^54^ 



where the Euchdean action is 

r/3 



^^[m] = / drf f £r[-i^A{m) ■ a^m] + E[m]|. (155) 

The massive m^ field is coupled to the massless field through the time derivative term and 
the constraint |m| = 1. For simplicity, we assume that the pseudospins are aligned along the 
X direction or we concentrate on a local patch of the spin texture, where the pseudospins are 
almost aligned along the x direction. Using the constraint, we can express m^ in terms of m'^ 
and m^. We only keep terms quadratic in m^ and m^. The monopole vector potential is given 
in a convenient gauge by A ^ (0, — m^, m^)/2. Previously we took the long wavelength limit 
to obtain the local action in real space since we were interested in results which would become 
exact in the limit of extremely smooth spin textures. Here we generalize our discussion in 
order to obtain an approximation for the full spectrum of collective modes. The calculations 
are identical to those detailed in Section |rV| except that we do not take the long wavelength 
limit. Here we report results only for the case of experimental interest, u = 1, where the 
expressions take a somewhat simpler form. For small rriz and rriy it is straightforward to 
integrate out the massive m^ field since the different momentum components decouple at 
the gaussian level 



e-5f,,[m-] 



+D,(g)|m^(a;„,g)|2 + D^(g)|m^(a;„,g)rt , (156) 



yielding 



^3 



Sffflm'] = E ( ^i^^^ + Dyi<l))\m^i-^n,q)\'. (157) 

We can read off the collective mode frequency from this expression: 

LUk = 87T^D,{k)Dy{k), (158) 

where, 



+ / 7|^^^(?) exp(-|g|V2) - / J^yA{q) exp(-|g|V2)e^'=^''}, (159) 



and 



Dy{k) = i- I ^VE{q) exp(-|g|V2)(l - e^'^'^) . (160) 

This spectrum agrees exactly with the dispersion relation in Refs. [ |],0, where the excita- 
tion energy for the single magnon state was obtained using the time-dependent Hartree-Fock 
approximation and using a single mode approximation combined with the Bogoliubov trans- 
formation. In the long wavelength limit the dispersion relation reduces to the result discussed 
above. 
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VII. KOSTERLITZ-THOULESS PHASE TRANSITION AND SPIN-CHANNEL 

SUPERFLUIDITY 

The linearly dispersing gapless mode discussed in the previous section and the absence of 
gapless charged excitations suggest that the system should show superfluid behavior in the 
pseudospin channel as has been noted previouslyufl. To make this suggestion more concrete 
we evaluate the linear response of the system to opposing electric fields in the two layers. 
Combining Eq. ||146|| with the continuity equation for the z component of pseudospin we find 
that 

(Tzz[q,u}) = — = . 161 

This conductivity gives the difference between oppositely directed charge currents which 
flow in the two layers in response to oppositely directed electric fields. (Note that no net 
current will be induced by such electric fields). In the long wavelength limit this leads to a 
frequency dependent conductivity equivalent to that of a superfluid with superfluid density 
proportional to pe'- 

a^^{uj) = — — 2— • (162) 

tun uj 

Note that by the Kronig-Kramers relation this conductivity must have a real part which 
is proportional to a. 6 function at zero frequency. We remark that the superfluid property 
requires not only a gap for charged excitations and the linearly dispersing collective mode but 
also a total oscillator strength which vanishes as g^ at long wavelengths. For the conductivity 
associated with response to electric fields in the same direction in the two layers, the total 
oscillator strength associated with intra-Landau-level excitations vanishesEil as q^ and the 
the collective mode has a gap at long wavelengths. These properties lead to a a{q, u) which 
vanishes as q^ in the long wavelength limit leading to the quantum Hall effect rather than 
to superfluidity. 

The above analysis is dependent on our analysis of the response functions of the double- 
layer system which does not include thermal fluctuations. At finite temperatures both 
meron and pseudospin collective mode thermal fluctuations have to be accounted for. As in 
other two-dimensional superfluids the linear response conductivity is still expected to vanish 
below a finite temperature. As we have discussed previously the low-energy excitations of the 
double-layer system consist of highly elliptical precessions of the spin about the direction of 
the local order parameter. It follows from a continuity equation and the equation of motion 
for the spin textures that the conserved number current corresponding to the z component 
of the pseudospin is related to the pseudospin by the usual minimal coupling prescription 
(see Appendix for further discussion of this result) 

Jzz = ^V<^ (163) 

n 

where the factor of two arises from the definition of J^z as the difference of the number 
currents in the two layers, and (f) is the projection of the pseudospin orientation onto the 
X — y plane, (p plays the same role as the phase of the superconducting order parameter 
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in a two-dimensional superfluid and pseudo-spin channel superfluidity will be coincident 
with pseudospin easy-plane ferromagnetisni. Likewise, as we discuss further below and in 
Appendix , the divergent superfluid conductivity implies zero pseudospin Hall resistivity 
below the KT temperature. 

In this system the Kosterlitz-Thouless (KT) phase transitions^ which separates the su- 
perfluid and normal states is expected to beE3 associated with the unbinding of meron pairs 
of opposite charge and opposite vorticity. In order to analyze the KT transition, it is better 
to work in real space, since the KT transition is controlled by the large-distance physics and 
the short- distance details can be effectively taken into account by the vortex core energy. In 
the limit of a large mass term {(3f' ^ pe) fluctuations in m^ will be small. After integrat- 
ing out the massive fleld and flnite frequency fluctuations in a gaussian approximation, we 
obtain an effective XY model 



S!ff = i:PE dh{V^)\ (164) 



'eff- 2' 

where (p is an angle denoting the direction of the spin in the XY plane. We know that this 
model undergoes a Kosterlitz-Thouless phase transition associated with the unbinding of 
bound vortex pairsEJ at approximately 

T, = ^PE. (165) 

For a 2D XY model corrections to this expression for Tc arise from flnite temperature spin- 
wave and vortex-antivortex polarization renormalizations of the effective spin stiffness at 
long distances. The magnitude of the corrections depends on details of the short- distance 
physics. For the 2D nearest-neighbor-coupling XY model on a square latticeES whose short- 
distance physics we believe to be similar to that of double-layer systems 

Te ~ 0.90pE. (166) 

(We note however that the double-layer system possesses charged excitations whose energies 
are larger but of the same order as the meron core energies.) In the present case Tc should 
be further reduced, especially as d/£ goes to zero, because of fluctuations out of the x — _v 
plane. However, numerical studies of the anisotropic 0(3) model on a square latticeE3 
demonstrate that these corrections are not very important except in the limit of extremely 
weak anisotropy. For example T^ is reduced by a less than a factor of two compared to 
Eq. ||166|| even for parameters which correspond to pE^'^/P ~ 30. Comparing with Fig.(^ 
we see that such weak anisotropics occur in double-layer systems only for d/£ < 0.3, a 
regime which is not experimentally accessible. (Quantum fluctuations only increase the 
anisotropy by decreasing pE-) In the following section we provide a quantitative estimate 
of the temperature scale expected for Tc by combining Eq. ||166|| with estimates of pe which 



include the effect of quantum fluctuations. As we discuss there, the principle source of 
uncertainty in the Tc estimate comes from attempting to estimate pE- 

It should be noted that the U(l) symmetry responsible for the XY model physics and the 
KT transition is robust under application of a bias voltage which induces a charge imbalance 
between the layers. Within the pseudospin picture, the spins simply tilt slightly out of the 
XY plane in the positive or negative z direction. This will reduce the component of the 
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spin in the XY plane and hence may lower the KT temperature slightly but it will not (for 
weak imbalance) destroy the ordering transition. This is in sharp contrast to the behavior 
for mmn states with m > n which require precise charge balance for their existence. This 



point is discussed in more detail in Section [VIII] below. 

In an ordinary superconducting film, the linear response conductivity is infinite below 
the KT transition temperature. More precisely the volt age- current relationship obeys: 

JPiT) 



V 



(167) 



where the exponent p(T) = 1 for T > Tkt and p(T) > 3 for T < Tkt- The exponent has a 
universal jump and changes discontinuously from 1 to 3 at the transition. The critical current 
is zero, i.e. there is always a finite voltage for any finite current. At any finite temperature 
Bogoliubov quasi-particles are present, however these do not produce dissipation because 
the current is carried by the superfluid and the voltage is (essentially) zero. The superfluid 
'shorts out' the normal fluid and so there is no electric field to produce motion of the normal 
fluid. The present problem is richer because of the existence of two layers in which the 
current can flow and because of the possible presence of a Hall voltage in the presence of 
currents. It is instructive to define phenomenologicalE3 transport coefficients (for the case 
of two identical layers) as follows: 
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Here the numerical subscripts label the two layers of the double-layer system. (Note that this 
phenomenology reflects the fact that the transport coefficients for the sum and differences 
of the currents or electric fields decouple.) 

In general there are four independent transport coefficients, allowing for contributions 
to both dissipative and Hall electric fields due to currents flowing in the same layer, p^^ and 



p^ , and the interlayer or dragp resistivities coming from currents flowing in the opposite 



layer, p^^ and p^y. The ratios of the electric field sums and differences to the current sums 
and differences are given by p^ + p^ and p^ — p^ respectively for both dissipative and Hall 
fields. The matrix of conductivity coefficients may be related to the matrix of resistivity 
coefficients by inverting Eq. ||168|| . We have argued above that the dissipative dc conductivity 
coefficient for current differences is infinite when T < Tkt- This implies that both the Hall 
and dissipative resistivities are zero for T < Tkt, «-e that pf^ = p^^ and p^y = p^y. Below 
the Kosterlitz-Thouless temperature there are only two independent dc linear transport 
coefficients. Moreover we know a great deal about these two transport coefficients because 
of the quantum Hall effect. When identical electric fields exist in each layer the total Hall 
conductance is nearly exactly quantized at low temperatures {a^y + a^y ~ ue^/h) and the 
dissipative conductance is activated {a^^ + a^^ ~ o"oexp(— A/2fcBr))- Here a^ is a non- 
universal constant and A is the gap for making charged particle-hole pairs; these charged 
objects are probably the meron pairs discussed previously. In terms of the linear resistivity 
matrix of Eq. ||168|| we conclude that for low temperatures 
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and 

pi = pI = ^- (170) 

Note that this imphes the occurrence of remarkable cross-talk phenomena. t^S For example 
if current is injected in the x direction in one layer but no current flows in the other layer, 
a quantized Hall field whose value corresponds to the density per layer will appear in both 
layers. 

Hol^ has recently considered the question of the stability of current flow in the spin 
channel. He finds that the existence of oppositely directed electric fields in each layer 
induces phase twists which can be relieved only by steady nucleation of topological defects 
in close analogy with textures induced in superfluid ^He. We agree that such effects will 
occur at finite current densities, however these results disagree with our findings above that 
the linear response a^x is infinite and p^y is zero. We believe that this discrepancy has 
two origins: the finite easy-plane anisotropy is essential to the analysis of this problem and 
the fact that there are subtleties associated with the question of the existence of uniform 
current carrying states in the lowest Landau level. [This latter point is discussed in detail in 
the Appendix.] Without easy plane anisotropy, the SU(2) symmetry prevents the XY order 
necessary for superfluidity. However since even with easy plane anisotropy, the XY order is 
only algebraic, the critical current density is zero as discussed above. For any finite current 
density there will indeed be dissipation due to generation of topological defects. 

VIII. EXACT DIAGONALIZATION STUDIES 

In this section we discuss some microscopically based investigations of the properties 
of double-layer systems at i^ = 1. We start by discussing some studies using the finite- 
size exact-diagonalization in the spherical geometryE^. As we have emphasized previously 
at d/i = 0, where intra- and inter-layer Coulomb potentials are the same, the electron- 
electron interaction term in the Hamiltonian is invariant under rotations in pseudospin space. 
Eigenstates of the Hamiltonian H can be simultaneous eigenstates of any one component 
of the total pseudo-spin operator. States with pseudo-spin quantum number S will have 
degeneracy 25* + 1. As discussed in Section |T| the ground state at i/ = 1 is a pseudo-spin 
eigenstate with S = N/2 where A^ is the number of electrons.) At finite layer separations, the 
inter-layer interactions will be weaker than intra-layer interactions and the broken symmetry 
in the ground state is reduced from SU(2) to U(l). Microscopically this corresponds to 
the fact that in the ground state, electrons within the same layer will be more strongly 
correlated than electrons in different layers in order to minimize the total Coulomb energy. 
For sufficiently widely separated layers we will provide evidence that the ground state no 
longer has broken symmetry. This problem has been studied before through examination of 
the ground state wavefunction,0E2l and a ground state level crossing was found to occur in 
the vicinity of d/i = 1.5. In this section we will discuss another attempt to estimate the layer 
separation at which this quantum phase transition takes place. In addition we will discuss 
rough estimates obtained for the charge-gap and for the parameters of the spin-texture 
energy functional. 
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The model we consider consists of two two-dimensional electron systems separated by 
a distance d. For the sake of definiteness the spread of the electron wave function in the 
perpendicular direction in each layer is neglected. (Such effects are easily accounted for if 
the geometry of a particular sample is known; we assume that it will normally be possible 
to define an effective layer separation for each sample.) We start by considering the case 
of greatest interest where each layer has the same number of electrons and define the zero 
of energy by placing neutralizing non-responding background charge backgrounds on each 
layer. We parameterize the intra- and inter-layer interactions in terms of the Haldane 
pseudopotential parameterg53 for this model. Calculations were performed for iV = 4, 6, 8, 10 
at a variety of layer separations. 

We first discuss finite size estimates of the chemical potential dependence at z/ = 1. 
For a finite system of electrons on the surface of a sphere the chemical potential jump is 
expected to occur when N = Norb = N^ + 1. Here Norb is the number of orbitals per Landau 
level when N^ flux quanta pass through the surface of the sphere. The chemical potential 
discontinuity is given by the limit d& N ^ oo oi 

A/i = E{N^, N + l) + E{N^, N-1)- 2E{N^, N), (171) 

where N^ = N — 1. In Fig. (pH]) , the finite size estimate of A/i for the u = 1 state is shown 



as a function of layer separation d for A^ = 8 and A^ = 10. Afi decreases continuously as d 
increases. The finite size corrections can be crudely inferred from comparisons of results for 
N = 8 and N = 10. For a system size increase from A^ = 8 to A^ = 10, A/i increases slightly 
for d < 1.5i, and decreases substantially for d > 1.5i. In all likelihood this result indicates 
that in the thermodynamic limit (A^ = oo), there is no chemical potential discontinuity 
for d >~ 1.5i. For layer separations larger than this, it seems likely that the properties 
of the double layer system should be similar to those of two isolated lavers with Landau 
level filling factor u = 1/2, which are of great interest in their own right.Ell The existence of 
an intermediate state, best viewed as a quantum disordered easy-plane ferromagnet, is an 
interesting speculative possibility. 

We have argued that the charged objects, at least away from d/i = consist of meron 
pairs with opposite vorticity. Since we expect the meron core energy and the pseudospin 
stiffness to both vanish when the order parameter vanishes, we expect that the chemical po- 
tential discontinuity should vanish at the critical layer separation where the zero-temperature 
phase transition to a disordered state occurs. Thus these results also provide an estimate 
of the layer separation at which the quantum phase-transition takes place. We can get an 
estimate of the finite-size errors in these results at small d/i = where the value of the 
chemical potential discontinuity is knownll3 exactly for u = 1. At z/ = 1 and d/i = 0, 
the ground state many-body wavefunction is a single Slater determinant so that there are 
no quantum fluctuation corrections to the Hartree-Fock value for the spin-stiffness given in 
Section |^. The chemical potential discontinuity is therefore given exactly for A^ = oo by 
A/i = (e^/£)(7r/8)^/^. This value is marked by an asterisk in Fig.(|Tl|) and compares well 
with the finite-size estimates. 

To get further insight into the system, we have attempted to estimate the dependence of 
the ground-state order parameter (M) and the easy-plane pseudospin magnetic susceptibility 
(x) on d/i for u = 1. We normalize the order parameter so that it has the value 1 at d/i 
=0. The order parameter is defined by 
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M^lim[ lim -il/N)^EiN^,N)\] (172) 

and the magnetic susceptibility is defined by 

X^lim[Jim -(l/iV)^E(iV^,iV)|] (173) 

where t is the inter-layer tunneling amplitude (which acts like a Zeeman field on the pseu- 
dospins). Some typicalO results for the dependence of the ground state energy on tunnehng 
amplitude are shown in Fig.(|T^ and Fig.(|T^. The solid lines in these figures are interpola- 
tions between energies calculated at a series of equally spaced t values not including t = 0. 
In a finite size system it is clear that the ground state energy always will have a quadratic 
dependence on tunneling amplitude until the total 'Zeeman' coupling becomes comparable 
to the smallest excitation energy of the system so that the order of limits in Eq. [|1 72|| and 



Eq. |[L73 ] is essential. The lowest energy excitations of the system will be long wavelength col- 



lective modes and we can expect to estimate the dependence of the zero-temperature order 
parameter and susceptibility on d/i from the dependence of the ground state energy on t in 
the regime where the t is small compared to e^/i but the Zeeman energy is large compared 
to the minimum excitation energy of the system. Using our result for the collective mode 
energy and the correspondence between linear momenta and angular momenta on a sphere 
{k ~ L/R where R is the sphere radius) we estimate that the minimum value of t we can 
consider is 



47rv/pg^ 

For u = 1, pe is strongly renormalized downward by quantum fiuctuations and we expect it 
to drop to zero as d/£ approaches ~ 1.5 while P starts from zero at d/i = 0, rises quadrati- 
cally and is less affected by quantum fluctuations. Comparing with Fig. (H) we estimate that 
for d/i ~ 1 and N = 10, tmin ~ I0~^{e^/i); tmin will approach zero for d/i — > 0. These rough 
estimates are consistent with the crossovers from quadratic to linear dependence on d seen 
in Fig.([T^) and Fig.(O). We have estimated the magnetization and susceptibility by inter- 



polating between calculated energy values and evaluating the derivatives at (i ~ 2.5 x 10 ^. 
The results are shown in Fig. (JT^ and Fig. ([T3|) . The magnetization decreases continuously 



with d and appears to vanish in the limit of large A^ at a critical layer separation consis- 
tent with the estimate obtained by looking at the charge gap. The flnite-size susceptibility 
estimate becomes large near where the order parameter vanishes as expectedO. 

The order parameter is an important parameter which will appear in the spin-texture 
energy functional in the presence of interlayer tunneling as discussed in Ref . [ |13[ . We discuss 



here attempts to estimate parameters (3 and pE- Fig.(|T^) shows the dependence on the 
layer separation of the increase in the ground state energy when one and two electrons are 
transferred from one layer to another. According to the spin-texture energy functional this 
energy should be given, in the limit of large systems, by 



AE = (3{d/i)i'^2Tx 



40 



(AiV)- 



(175) 



The finite size exact diagonalization calculations were performed for d/i = 0.5, d/C. = 1.0, 
and d/i = 1.5 and the values of /?£^/(e^/£) inferred by comparing with Eq. ||175|| were 5.1 x 
10~^, 1.7 X 10"^ and 4.0 x 10~^ respectively. These values compared to 2.0 x 10~^^, 4.0 x 10~^ 
and 6.0 x 10~^ respectively for the Hartree approximation and 7.3 x 10~^, 2.2 x 10~^ and 3.8 x 
10^^ for the Hartree-Fock approximation results derived earlier. The quadratic dependence 
derived in the Hartree-Fock approximation is apparent in the exact diagonalization results. 
The exact diagonalization estimates for (3 demonstrate that quantum fluctuation corrections 
to this quantity are quite small and that the Hartree-Fock results shown in Fig.(^ are quite 
accurate. An important aspect associated with the broken symmetry in the ground state 
at z/ = 1 is the fact that the chemical potential discontinuity is not strongly influenced by 
the transfer of charge from one layer to the other. In the presence of a bias potential the 
ground state pseudospin is tilted out of the x — y plane but the system still has a broken 
U{1) symmetry. This situation contrasts with the case of the double-layer quantum Hall 
effect which occurs at total Landau level filling factor z/^^ = 1/2 where the chemical potential 
discontinuity occurs only near equal layer densities.ll3 

In Fig. (p!7|) and Fig.(|T8|) we show the dependence of the lowest excitation energy on 
wavevector for a finite system of electrons on a torus with periodic boundary conditiongiHl 
for d/l = 0.5 and d/i = 1.0 respectively. At d/i = 0.5 clear indications of the linearly 
dispersing collective mode are evident. The velocity of the linearly dispersing collective 
mode is related to the parameters of the spin-texture energy by 

Ait 



E(k) = ki—J2(3pE. (176) 

Reading off the velocity from Fig. (|17D and using the value of /3 obtained above from the 
exact diagonalization calculation we estimate that for d/i = 0.5, pe ~ 1.5 x 10~^(e^/£) 
compared to the Hartree-Fock value pE = 1.22 x 10~^. For d/i = 1.0 the linear dispersion of 
the collective mode is less evident and it is more difficult to estimate the velocity accurately. 
We have chosen to estimate the velocity from the energy at the smallest finite-size system 
wavevector as indicated in Fig.(|T8D from which we obtain pE ~ 3.3 x 10~^(e^/£) compared 
to the Hartree-Fock value pE = 4.22 x 10"'^. At this value it is evident that quantum 
fiuctuations are decreasing the value of pe as the critical layer separation is approached. 

To obtain a rough quantitative estimate of the layer separation dependence of the 
Kosterlitz-Thouless transition temperature, we calculate the leading order corrections to 
the Hartree-Fock ground state energy as a function of parallel magnetic field B\\. The renor- 
malized spin stiffness p^ is related to the dependence of the ground state energy on field by 
the following equations 

p| = ^E^(g)|Q^o (177) 

where Ea is ground state energy per area and Q = dB\\/i'^Bp. Given p^ and the fact 
that the mass term is not strongly infiuenced by quantum fiuctuations, we can estimate the 
effect of fluctuations in ttt,^ on the transition temperature by using results from numerical 
simulationaB of the classical 0(3) 'easy-plane' non-linear sigma model. We flndE^ the de- 
pendence of Tkt on d/i which is illustrated in Fig.( |T^. As the layer separation approaches 
zero, U{1) symmetry is enhanced to 0(3), and fluctuations of m^ eventually become suf- 
flciently important that Tkt decreases and approaches zero. On the other hand for layer 
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separation close to the critical one, phase coherence between the layers is destroyed only by 
quantum fluctuations. Our results for T^t are in good agreement with the d/i dependence 
of the temperatures at which features in the dissipative resistance of double-layer systems 
have been seenEZl by Lay et. al.. These features might be associated with the vanishing 
charge gap expected to occur in parallel with the KT transition. 

IX. CHERN-SIMONS-LANDAU-GINZBURG THEORY 

An useful alternative way of understanding the physics of the quantized Hall effect is 
based on the concept of composite bosonalaEl. Much of the physics of the single layer quan- 
tum Hall effect can be described in terms of the Chern-Simons-Landau-Ginzburg theoryE^El. 
One starts with the observation that the problem at odd denominator filling fraction can be 
mapped exactly to a problem of bosons with an odd number of statistical flux tubes attached 
to them. At the mean field level, the statistical flux cancels the external magnetic field and 
one obtains a boson superfiuid. Treating the fluctuations above this mean field within 
the RPA approximation restores the gap and one obtains an incompressible bosonic liquid. 
Both Laughlin's wave functioned and the long wavelength algebraic off- diagonal-long-order 
correlation function^ can be derived explicitly from the CSLG theoryEll. 

The CSLG formulation of the single layer quantum Hall system has been extended both 
to the case involving electron spinll^O and to the double layer caset^'OflS by a number 
of authors. The Chern-Simons theory of the double layer system based on the fermionic 
representation has also been constructedEHl. In the case of double layer systems, there is 
both a statistical gauge interaction of the composite bosons within the layers and between 
the layers. The action for this problem is given by 



C^ = (p\{idt + A° - a°)0^ - ^— ^ (^_ V + A - ma^ - na^) </ 

'£y5p^{x)VAix-y)Sp^{y) + {]^l)- [ £y Sp^{x) Ve{x - y) Sp^{y), (178) 



2 

where (T~^i) indicates the first four corresponding terms with up and down labels inter- 
changed. Here m is an odd integer. In the absence of tunneling, the particles in the two 
different layers are distinguishable, the relative phase winding between them can be either 
or IT, therefore, the integer n can be either even or odd. At the mean field level, the 
equations of motion are given by 

V X a^{x) = 2iTp^{x) ; A = ma.^ + naj^. (179) 

When the electron densities of both layers are equal, we see that the the filling factor has 
to be z/ = 2/(171 + n) for these equations to be satisfied. 

From this CSLG formulation Halperin's wave functional for the [mmn) state can be 
derived in a fashion similar to the single layer caseEll. One can decompose the complex 
boson field in terms of the amplitude and the phase part. 



(t)M = ^P + 5p.e''"^^\ (180) 
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The Chern- Simons gauge field induces a long-ranged logarithmic density-density interaction, 
giving rise to the following effective Hamiltonian at the quadratic level: 



2 m* q 



— j ((m^ + n2)7r,(g) j^^^-q) + Ammi^{q) T^i{-q))+q^ 9^{q) 9^{-q) 



(181) 

where vr^ is proportional to the density 6pa^{x) and is the conjugate variable of the phase, 
i.e. 

[Oaiq), vr.'(g')] = -^ S<^cr'S{q + q'). (182) 

This Hamiltonian is a direct sum of harmonic oscillator Hamiltonians for each wave vector 
q. By forming the sum and the differences of these operators, one can easily diagonalize 
the Hamiltonian and find the ground state wave function in terms of its dependence on the 
generalized coordinates (densities) vr^: 

r n f 1 .5-^ 27r(?7i + n) , , 2-11 (m — n) 1 

One can express the density operators in terms of the ordinary first quantized coordinates 
of the particles 

< = ^ E {e^'^' + e''^'^ - 2p} ; ''^ = ^^ {^^''" - ^^''" } ' (^^4) 

and obtain the first quantized wave function 

^o(a;J,a;]) = H l4 - x]r \xl - xJT \xl - xjre-^^-.-l"^''. (185) 

i<j 

This is nothing but Halperin's wave function for the {mmn) state expressed in terms of the 
composite boson variables. One can easily perform a singular gauge transformation back to 
the original electron coordinates and find explicitly Halperin's wave function expressed in 
terms of the original electrons. 

From the effective Hamiltonian ( |181| ) one can also derive all static correlation functions in 
the long wavelength limit. One example is the analogue of the ODLRO correlation function 
first introduced by Girvin and MacDonaldEs for the single layer case. From ( |181| ) one easily 
obtains 

< e^^+(^)-^^+(s^) >= \x - y|-('"+") ; < e*^-(^)-^^-(s^) >= \x - y\-^"'-''\ (186) 

Similarly, one can obtain the static density correlation function in the long wavelength limit 
(for m ^ n) 

q^ o^ 

< vr+ 7r+ >= — ; < 7r_" ttZ^ >= — -. (187) 

'' '^ 2Tx{m + n) '^ '^ 27r(m-n) ^ ' 
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We therefore conclude that the CSLG theory correctly gives the static properties of the 
double layer systems in the long wavelength limit, and they agree extremely well with the 
results of the microscopic calculations based on Halperin's wave function. At the level of the 
static correlation functions, the degree of agreement is similar to the single layer caseLjH. 
However, there is some discrepancy in the collective mode spectrumEl between the CSLG 
theory and the microscopic theory obtained using the projected single mode approximation. 
Within the CSLG theory, the collective mode spectrum can be obtained by studying the 
Gaussian fluctuations of the bose order parameter, and one obtains the following spectrum 
for the in-phase and the out-of-phase collective modes 

in — n 

UJ+ = nujc ; uj- = riuJc , (188) 

m + n 

in the long wavelength limit. One sees that the in phase mode agrees exactly with the 
prediction of the Kohn's theorem, since the in-phase magnetic translation is a good symme- 
try. The out-of-phase mode agrees exactly with the SMA spectrum obtained using the full 
density operator pq 

<P^'^''^P-^'^> = nu:^^^. (189) 

< Pqp-q > m + n 

This result is not surprising, since this SMA formula involves only the static correlation 
functions of the full density operator, and these are obtained correctly within the CSLG 
theory, as shown above. The drawback of the CSLG theory lies in the fact that there is 
no sensible way of obtaining the projected density operators. In the limit of large Landau 
level spacing, the SMA using the full density operator is not adequate, and the projected 
operators must be used. Explicitly projecting onto transitions from the iV = to the N = 1 
Landau levels, one obtains^ 

uj^ = hu, ;co_ = hco,-l^^q\''h^iq). (190) 

for the inter-Landau level transition modes. From a microscopic point of view, this dis- 
crepancy should be resolved. However, one can also take a phenomenological point of view, 
and fix the mode frequencies obtained within the CSLG theory as parameters fitted to the 
projected SMA calculations. 

It is clear for the case m = n that the above formulation breaks down. The true wave 
function is not simply given by Eq. ( |183| ) which lacks all correlations in the 7r~ channel. 



Instead one must build in Jastrow correlation of the form discussed in Section ^ in order 
to obtain the correct linear dispersion of the long wavelength spin fluctuations. 

For the special case of the {mmm) states, the SU{2) symmetry of the problem has been 
exploited by Lee and Kaneli^, and they formulated a slightly different version of the Chern- 
Simons theory that is very useful for uncovering the spin charge connection. Their original 
motivation was to understand the spin unpolarized quantum Hall effect. However, with a 
simple change from spin to pseudospin, one can easily apply this formalism to double-layer 
systems (in the zero-separation limit where the SU(2) symmetry is preserved). In the case 
of the {mmm) states, one only needs to introduce a single gauge field, and the action is 
given by 
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L^ = (Plitdt + A° - a°)0. 



2m* 



1 

- V + A — ma 



Air 



~ 2J ^^^ "^Palx) Vo{x - y) 6p^{y). 
In this case, the mean field equations are given by 

V X a(x) = 27rp^(x) ; A 



m,a 



(191) 



(192) 



and one sees easily that these equations are satisfied when the filling fraction is z/ = 1/m. 
We now decompose the bose fields in the form 



p + 6p^ (j) z„ 



^cr ^cr -L • 



(193) 



Here is a complex number of unit magnitude and z^j represents the spinor variable and is 
related to the unit vector order parameter m definedjjreviously by m 



^fl ^ ^v ^v 



Performing the standard duality transformationE^, one obtains the following effective 
action in the dual representation 



m, 



£ = 27r 6„ {Jl + r^) + — (9o b^f 

2p 

27r 



Air m 
- I (e"^ d^ 5hp) V{x - i/)(£"'^' do., 5hp. 



27r{J^ + J^) e""^ a„ 6b 



In \x — y\ 



MJo + Jo] 



27r 



£°'^' 9„/ 6bs, 



where Sba = ba + ^^a and the gauge field 6^ is defined by 



Ju. 



^fivp^u^p 



with J a being the three current of the fiuid, 



1 



Jl = 7r-^pupdu{(l)dp(j)) 



2m 



(194) 



(195) 



(196) 



is the vortex three current and 



^ 2m 



(197) 



is the Skyrmion three current. Note that the word vortex in this context refers to vortices 
in the bosonic Chern-Simons field, and should not be confused with vortices in the spin field 
discussed in the other sections in terms of the pseudospin XY ferromagnet analogy. 

From this dual action, several important results follow. First of all, one sees that there 
is a long ranged logarithmic interaction of the topological density with itself Ptop(^) = 
Jq + Jq — i^"^ da Sbj3. In the low energy sector, the only excitations which can be created 
are those which have no net topological charge, i.e. J d'^rptop{r) = 0. Therefore, these 
elementary excitations can be classified into three categories. A vortex excitation has 



JV 



-e"'^ da 6bp. 



(198) 



45 



Therefore, it carries charge ±i>e depending on the sign of the vorticity. A Skyrmion excita- 
tion has 

J^ = -e^^ d^ 5hp (199) 

and it also has charge ±z/e, depending on the sign of the Skyrmion number. This is the 
equation in the Chern-Simons theory which relates the spin and charge, which was first 
noted by Sondhi et alt3. We see that this relation is exactly the same as the one obtained 
in Eq.(l^) from a microscopic calculation. One can also form a bound state between these 
objects so that the net charge is zero. 

r, = -r, (200) 

In this case, the vorticity cancels the Skyrmion number exactly, these objects are neutral. 

The statistics of various spin textures can also be worked out explicitly from the dual 
action, following an approach used by Lee and Zhangta in the case of the CSLG theory for 
the single-layer QHE. From Eqs.( |194| ), one sees that the first term 2iiba{Ja + •Jo) couples the 
Skyrmion density Jg to the dual gauge filed ba- The coupling to the averaged flux (ha) gives 
the dynamics of the spin degrees of freedom discussed below. There is also a coupling to the 
fluctuating part, 56/3, which is given by Eq.( |199| ). Therefore, this gives rise to a statistical 
interaction between the Skyrmion density and the Skyrmion current of the usual form 

pj d'xd'yJ^e'^^^^^^r^iy). (201) 



Another important consequence of the dual action in Eq.( p.94| ) is the form of the effective 
spin action. In terms of the CP^ fields, it is given by 

K 

C = ip %dtz„) - — (|9z<^|^ + (z^9z<^)^) - iJ^(z^9^z<^). (202) 

Without the last term which couples the spin and charge degrees of freedom, this action 
can be transformed into the angular variables m. = z a z, and its form agrees exactly with 
the effective spin action in Eq.(^) derived earlier from the microscopic calculation. The 
microscopic calculation was carried out with the assumption that the charge degrees of free- 
dom is massive. Under this assumption, integrating out the charge degrees of freedom in 
Eq. ( P02| ) will only produce a long ranged Coulomb interaction between the topological den- 



sity. However, strictly speaking, this assumption is true only at zero temperature. Eq. ( |202| ) 
is more generally valid even at finite temperature where the charge degrees of freedom is 
gapless. The effect of the gapless charge degrees of freedom on the spin dynamics and the 
Kosterlitz-Thouless transition is still to be explored. At this level, the coefficient K = p/m* 
derived from the Chern-Simons theory depends on the mass of the electron, rather than 
the Coulomb interaction as it should. This is a general feature encountered in all Chern- 
Simons theories. One can view this coefficient as a parameter and argue that higher-order 
corrections will bring it into agreement with microscopic theories. The coefficient of the time 
dependent term is independent of the mass, and agrees exactly with the result of microscopic 
calculations. At zero temperature, the charge degrees of freedom have a gap, integrating 

46 



them out would only give rise to a higher derivative coupling between the spin variables. 
Therefore, at least for zero temperature, the effective spin action given here is sufficient. 

Finally we note that, as the layer separation exceeds the critical value d*, quantum 
disordering will cause merons to proliferate. The path integral configurations for the meron 
world lines are very similar to those of the vortices which disorder the 3D XY model. However 
the universality class of the transition is different because of the coupling to the Chern- 
Simons field.E3 



X. SUMMARY 

We have presented here a theory of the spontaneous development of interlayer phase 
coherence in double-layer quantum Hall systems at various filling factors. Using a pseudospin 
language we have shown that the system is equivalent to an easy plane itinerant ferromagnet 
with an unusual spin-charge connection. There is a zero-temperature phase transition to 
a quantum disordered phase if the layer separation exceeds a critical value d > d*. For 
< d < d*, the system is predicted to exhibit a finite temperature Kosterlitz-Thouless 
transition, in the absence of interlayer tunneling. Our theory is expected to apply to any 
filling factor at which there is an incompressible state which is not a pseudospin singlet. 
Here however we have concentrated primarily on the case of filling factor z/ = 1. 

In a companion papertJ we will discuss the new phase transitions that occur in the 
presence of tunneling and parallel magnetic field. We will also make contact with the recent 
experiments of Murphy et al.Ej which appear to have observed one of these phase transitions. 
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XII. APPENDIX: GAUGE INVARIANCE AND UNIFORM CURRENTS IN 

LANDAU LEVELS 

The (static) effective action we have derived for the case of easy plane anisotropy has as 
its leading gradient term 

S = ^PeIV^I'. (203) 

The phase (p describes the local spinor orientation 

[ e+../2 j ^ (204) 

and the 'charge' conjugate to ip{r) is >S'^(r), which gives the local (physical) charge difference 
between the layers. In order to study the gauge symmetry in this problem, it is convenient 
to introduce charge and pseudospin gauge fields 

A± = A| ± A|, (205) 

where A| and A^ are the electromagnetic vector potentials in each of the layers. The order 
parameter ip is gauge neutral with respect to A+ since it corresponds to the condensation 
of a (physical) charge-neutral operator 

e^^Woc(^|(r)^i(r)). (206) 



However a gauge change 



modifies the wave function 



A_ — A_ + 2^Vx- (207) 



M_e.p(a,_S') (:*;), (208, 



which means that the action has the usual minimal coupling form 

1 ,_ e 



5 = ^P£;|V</^ + ^A„|^ (209) 

and the pseudospin current is given (for A_ = 0) by 

Jzz = ^V0 (210) 

which is identical to the result derived in Eq.( |163| ) using the equation of motion for the 
projected density. As mentioned previously in Eq.( |15U| ), the superfluid mode has oscillator 
strength proportional to q^ as expected for an ordinary superfluid. However the coefficient 
is proportional to pe and hence non-universal. 
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The minimal coupling argument given above is quite correct, however there are remark- 
ably confusing subtleties lurking just beneath the surface. For V0 equal to a constant, 
Eq.( piO| ) implies the existence of a uniform zero wavevector current in the LLL. This is 
paradoxical as we now discuss. 

The question of the form of the current operator in the LLL is a subtle one which has been 
considered by several author s.IHI1eJE2 The difficulty lies in the fact that the true unprojected 
current operator has no matrix elements within the LLL. It is purely off-diagonal, taking for 
example, the A^ = Landau level into the A^ = 1 Landau level. Hence it would appear to be 
impossible to have low energy currents at zero wave vector, despite our previous derivation 
using the fact that the projected densities do not commute. The resolution of this paradox 
can be seen in the following simple example. 

Consider the SU(2) invariant z/ = 1 case and let the ground state be fully polarized up. 
Restricting H to the LLL, this state is an exact eigenstate. The exact single-magnon excited 
states 

^k = ^^, (211) 

are labeled by a conserved momentumc9 due to the fact that they are gauge neutral (with 
respect to A+. They correspond to spin flip magneto-excitons and have a velocity 

V = -VkCk X I = -Vk^(^^k) X z, (212) 

Ti h 

(where ^(ri — r2) is the particle interaction) and hence would seem to have a finite pseu- 
dospin current at g = u; = 0. If we use the true unprojected current operator, then we 
must take care to include LL mixing in \&k- First order perturbation theory shows that such 
mixing is small and proportional to ilWl/huJc- However the current 

has matrix elements connecting adjacent Landau levels proportional to iuc, so that Landau 
level mixing, though small, gives a crucial contribution to the current, independent of the 
smallness of the mass m. Carrying out the perturbation theory in detail, yields results in 
complete agreement (to first order in ilWl/huc) with the minimal coupling considerations 
and the equations of motion methods discussed above. 
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FIGURES 

FIG. 1. Schematic conduction band edge profile for a double-layer two-dimensional electron 
gas system. 

FIG. 2. Illustration of magnetic translations and phase factors. When an electron travels 
around a paralellogram (generated by TqTkT-qT-k ) it picks up a phase (j) = 27r^ = q^ where <I> is 
the flux enclosed in the paralellogram and $o is the magnetic flux quantum. 

FIG. 3. Anisotropy mass (in unit of e'^/ei^) and easy-plane spin stiffness (in unit of e'^ /ei) as 
a function of layer separation for u = 1. These results do not include quantum fluctuations which 
are important at finite layer separation. 

FIG. 4. Anisotropy mass and easy-plane spin stiffness as a function of layer separation for 
ly = 1/3. These results do not include the dependence of quantum fluctuations on the layer 
separation. 



FIG. 5. Illustration of a skyrmion on a sphere. 

FIG. 6. Illustration of localized quasiparticle and quasihole excitations at v = 1. 

FIG. 7. The number of spins flipped in the ground state versus the system size (Landau 
level degeneracy) when a single electron is added to the u = 1 incompressible state on a torus. 
These results demonstrate that the ground state contains a single skyrmion spin-texture whose size 
is determined by the competition between minimizing the frustration required by the boundary 
condition and the Coulomb energy. 

FIG. 8. Microscopic skyrmion energy vs. the scale size A. The trial wavefunction interpolates 
continuously between a spin-polarized quasihole at A = and a smooth skyrmion spin-texture for 
A —> oo. 

FIG. 9. Illustration of merons (vortices). The spin configuration in the core region tips 
smoothly out of the XY plane making this object essentially one-half of a skyrmion. 

FIG. 10. Illustration of a meron pair with opposite vorticity and like charge. We propose that 
under appropriate circumstances these objects can form the lowest energy charged excitations in 
the system. 

FIG. 11. Chemical potential discontinuity A^u as a function of layer separation d for v = \. 
The results are for system sizes of eight and ten electrons. The '*' mark shows the value of A/Li at 
d = and N = oo according to the skyrmion theory. 
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FIG. 12. Ground state energy versus tunneling amplitude at small tunneling amplitudes at 

dji = 0.2 for iV = 4, iV = 6 and iV = 10. 

FIG. 13. Ground state energy versus tunneling amplitude at small tunneling amplitude for 
AT = 10 at d/i = 0, d/l = 0.2 and d/i = 0.6. 

FIG. 14. Finite size estimate for the ground-state order parameter as a function of layer sepa- 
ration for V = 1. The results are for system sizes of six, eight, and ten electrons. 

FIG. 15. Finite size estimate for the ground-state magnetic susceptibility as a function of layer 
separation for v = 1. The results are for system sizes of six, eight, and ten electrons. 

FIG. 16. Dependence of the ground state energy level on d/i for one and two electrons 
transferred between wells. 

FIG. 17. Wavevector dependence of low energy excited states for i/ = 1 , d/i = 0.5 and 
N = 10. These results can be used to estimate the quantum renormalized spin-stiffness. 

FIG. 18. Wavevector dependence of low energy excited states for u = 1 , d/i = 1.0 and 
N = 10. These results can be used to estimate the quantum renormalized spin-stiffness. 

FIG. 19. Estimate of Kosterlitz-Thouless transition temperature for B = 4.364 Tesla so that 
e^ jikB ~ 106 Kelvin. These estimates include quantum renormalizations of the spin-stiffness and 
corrections to the XY model due to finite anisotropy strength. 
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